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ABSTRACT 


The  results  of  an  experimental  Investigation  of  airfoil  dynamic 
stall  involving  large  amplitude  pitching  motions  are  described.  The 
discussion  is  focused  on  constant  pitch  rate  motion  histories.  Measure¬ 
ments  of  unsteady  surface  pressure  distributions  for  both  NACA  0012  and 
NACA  64jA012(13)  profiles  have  been  employed  to  infer  the  onset  and 
evolution  of  an  energetic  leading  edge  separation  vortex  over  a  wide 
range  of  flow  and  motion  conditions.  These  data  have  also  been  inte¬ 
grated  to  provide  estimates  of  the  time  histories  of  the  lift,  pressure 
drag  and  moment  coefficients.  The  effectiveness  of  the  various  motions 
for  achieving  lift  enhancement  has  been  determined  through  the  introduc¬ 
tion  of  a  dimensionless  Impulse  parameter  which  accounts  for  both  the 
magnitude  and  duration  of  the  additional  lift  increment  during  the  post¬ 
stall  period. 


I .  INTRODUCTION 


The  problem  of  airfoil  dynamic  stall  has  been  researched  extensively 
since  it  became  a  prominent  topic  during  the  late  1960's.  Most  studies 
of  this  interesting  aerodynamic  phenomenon  have  been  motivated  by  the 
desire  to  understand  the  source  and  behavior  of  torsional  helicopter 
blade  oscillations  and  their  effect  on  high  speed  rotor  performance.  As 
a  consequence,  most  of  the  available  results  have  been  generated  for  a 
range  of  conditions  and  lifting  configurations  which  attempt  to  emulate 

the  rotor  flight  environment . ^"9 

The  present  investigation  is  concerned  with  the  evolution  of  the 
dynamic  stall  flowfield  from  a  more  fundamental  perspective  and  with  a 
decidedly  different  motivation.  Specifically,  it  is  focused  on  several 
aspects  of  the  process  including  unsteady  vortex  lift  generation  with  an 
orientation  toward  potential  productive  applications. 

The  essential  features  of  the  dynamic  stall  flowfield  created  by 
harmonic  airfoil  oscillations  in  an  otherwise  steady  and  uniform  free- 
stream  have  been  described  by  McCroskey  and  his  coworkers . ^“3  The 
behavior  of  the  flow  for  what  has  been  termed  the  "deep  stall"  case  is 
of  special  significance  here.  For  large  amplitude  oscilletions  in  pitch, 
the  flow  near  the  airfoil  surface  is  observed  to  remain  "attached"  during 
the  airfoil  upstroke  to  a  point  well  beyond  the  static  stall  incidence 
angle.  This  delay  in  incipient  separation  is  followed  by  the  formation 
of  an  energetic  vortex  structure  near  the  airfoil  leading  edge  which 
grows  with  time  and  convects  downstream  as  the  airfoil  motion  proceeds. 

A  steep  suction  peak  in  the  instantaneous  pressure  distribution  on  the 
upper  surface  is  observed  to  remain  coincident  with  the  location  of  the 
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vortex.  The  resulting  high  lift  transient  and  subsequent  severe  nose 
down  pitching  moment  (sometimes  referred  to  as  "moment  stall")  are  now 
familiar  to  many  aerodynamists . 

Several,  more  general  observations  can  be  distilled  from  an  examina¬ 
tion  of  these  results.  The  first  of  these  is  that  this  class  of  flow 
behavior  must  be  viewed  and  treated  as  a  "process;"  that  is,  these  unsteady 
flows  are  comprised  of  a  sequential  time  history  of  events  whose  character¬ 
istics,  at  any  instant,  are  the  result  of  a  complex  interaction  between 
the  external  flow  environment,  the  past  flow  history  and  the  continuously 
changing,  time  dependent  boundary  conditions.  More  profoundly,  these 
results  suggest  that  the  term  "flow  separation"  be  redefined  to  accommo¬ 
date  the  notion  of  flowfield  reorganization,  if  only  on  a  transient 
basis,  which  is  apparent  in  many  of  these  unsteady  flow  problems.  In 
this  context,  the  term  separation  is  not  meant  to  infer  a  spatially 
local  event, but  rather  to  describe  the  global  flow  behavior  associated 
with  a  particular  configuration  and  its  time  history  of  motion. I**1 

In  this  light,  the  unsteady  separation  process  associated  with 
dynamic  airfoil  excursions  might  be  viewed  as  a  possible  vehicle  for 
extracting  freestream  kinetic  energy  for  the  purpose  of  enhanced  aero¬ 
dynamic  lift  generation.  Of  course,  this  process  must  result  in  an 
increase  in  the  overall  airfoil  drag  coefficient,  but  there  may  be  appli¬ 
cations  where  this  concomitent  high  lift,  high  drag  phenomenon  could  be 
productively  exploited,  e.g.,  for  flight  vehicle  maneuver  enhancement. 

Special  emphasis  is  placed  on  the  role  of  the  dynamic  separation 
vortex  in  the  unsteady  lift  generation  process.  In  this  regard,  it  is 
of  interest  to  consider  those  conditions  which  might  result  in  increased 
energization  of  the  vortex  and  the  separate,  but  related,  problem  of 
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optimizing  the  entrapment  time  of  the  vortex  structure  near  the  airfoil 
surface.  The  concept  of  employing  an  externally  "trapped"  vortex  for 
enhanced  lift  is  not  original  here, 12  but  the  idea  of  generating  and 
controlling  the  characteristics  of  an  unsteady  separation  vortex  for 
this  purpose  is  unique.  (This  idea  is  also  being  explored  independently 
by  Chow  and  his  coworkers. 13)  The  problem  can  be  viewed  as  one  of 
optimizing  the  time  dependent  boundary  conditions  to  maximize  the  lift 
augmentation. 

As  a  first  step,  experiments  have  been  devised  to  provide  new  insight 
into  the  influence  of  the  motion  time  history  on  the  dynamic  stall  flow- 
field  and  resultant  airfoil  loading.  This  report  describes  and  discusses 
the  results  of  a  preliminary  investigation  which  focuses  on  anharmonic 
unsteady  airfoil  motions  with  large  amplitude  variations.  Instantaneous 
pressure  measurements  have  been  obtained  over  upper  and  lower  airfoil 
surfaces  and  used  in  the  reconstruction  of  phase-averaged  pressure 
distributions.  These,  in  turn,  have  been  integrated  spatially  to  estimate 
force  and  pitching  moment  time  histories.  All  measured  time-dependent 
variables  have  been  correlated  with  the  appropriate  airfoil  motion  history 

The  study  was  restricted  to  conventional  (i.e.,  rigid  shape)  airfoil 
geometries  due  to  limitations  in  available  experimental  capabilities. 

In  principal,  however,  it  could  have  been  extended  to  include  the  effects 
of  arbitrary,  time  dependent  airfoil  shapes.  The  concept  of  employing  a 
"flexible"  lifting  surface  in  combination  with  equivalent  rigid  body 
motions  to  optimize  unsteady  separation  effects  is  an  intriguing  one. 

Although  acknowledged  to  have  a  potentially  profound  effect  on  the 
dynamic  stall  process,  the  "motion  time  history"  is  one  of  the  least 
understood  contributors  to  the  overall  flow  development . ^ >8  Several 
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prior  studies  have  been  performed  in  an  attempt  to  assess  the  signifi¬ 
cance  of  this  facet  of  behavior.  Besides  pure  pitching  motion,  unsteady 
free  stream  (or,  surge)  oscillations  have  been  investigated  for  fixed 
airfoil  incidence  conditions . ^ , 15  Limited  studies  of  combined  pitching 
and  surging  oscillations  have  also  been  conducted. Plunging  airfoil 
motion  has  received  recent  attention. 18  Despite  these  efforts,  a  unified 
and  physically  consistent  picture  of  the  effects  of  motion  time  history 
has  yet  to  emerge.  This  is  due,  in  part,  to  the  extremely  Large  range 
of  parametric  variables  which  must  be  considered  to  adequately  characterize 
these  effects. 

Some  progress  has  been  made  for  the  case  of  pure  pitching  oscilla¬ 
tions  about  a  point  interior  to  the  airfoil.  Early  experiments  involving 
enharmonic  pitch  motions  revealed  the  importance  of  pitch  rate  (a)  as  a 
dominant  variable.  This  was  suggested  by  Kramer^  as  early  as  1932  and 
confirmed  in  relatively  recent  measurements  by  Ham  and  Garelick^  and 
Kawashima,  et  al.20 

The  present  experimental  investigation  was  designed  to  complement 
prior  studies  by  extending  the  range  of  conditions  of  known  important 
parameters  (such  as  pitch  rate)  and  examining  the  effects  of  motion 
variables  and  combinations  of  variables  which  have  not  yet  been  adequately 
studied.  A  unique  model  oscillation  system  described  in  the  next  section 
provided  a  capability  to  study  motion  in  two  degrees  of  freedom  simul¬ 
taneously  -  one  in  rotation  (pitch),  and  one  in  translation  (either 
surging  or  plunging  motion).  To  date,  experiments  have  been  performed 
to  examine : 

(a)  constant  pitch  rate  motion  (a  =  constant) 

(b)  constant  pitch  rate  motion  followed  by  sudden  reversal  in 
pitch  direction 


(c)  constant  pitch  acceleration  (<*  -  constant) 

(d)  harmonic  pitch  oscillations  with  variable  termination  angle 

(e)  combined  constant  pitch  rate  and  constant  axial  acceleration 

(f)  combined  constant  pitch  rate  and  constant  axial  deceleration, 

and 

(g)  combined  constant  pitch  rate  and  constant  axial  translation 
rate  (impulsive  surge). 

Because  of  the  large  parameter  space,  the  amount  of  data  involved 
in  assessing  even  this  limited  class  of  motions  over  a  reasonable  range 
is  quite  extensive.  This  report  is  therefore  confined  to  a  discussion 
of  the  measurements  of  constant  pitch  rate  only.  Because  of  the  impor¬ 
tance  of  this  variable  in  establishing  the  character  of  the  dynamic 
stall  flowfield,  a  disproportionately  large  number  of  measurements  were 
obtained  for  this  condition  [cases  (a)  and  (b)].  The  effects  of  other 
airfoil  motion  combinations  will  be  discussed  in  a  later  report. 

The  parametric  effects  described  in  this  report  include  the  influence 
of  nondimens ional  pitch  rate  (described  through  the  variable,  k  = 
and  the  pitch  amplitude , a max.  The  pitch  rate  range  was  varied  over 
nearly  three  orders  of  magnitude  (0.001 <  k  <  0.2).  All  motion  was 
initiated  at  an  initial  airfoil  incidence  angle  of  zero  degrees.  The 
amplitude  variation  extended  to  extremely  high  incidence  angles 
(°max  =  60°)  well  in  excess  of  values  employed  in  any  known  prior 
studies. 

The  influence  of  airfoil  shape  was  also  investigated  to  the  extent 
that  two  different  airfoil  models  were  employed  in  the  study.  A  6-inch 
chord  (c  =  15.24  cm)  NACA  0012  airfoil  model  was  used  in  an  initial  test 
sequence  and  later  complemented  by  an  8-inch  chord  (c  =  20.32  cm)  NACA 
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64]A012(13)  airfoil  model.  The  use  of  these  two  different  shapes  provided 
a  limited  measure  of  the  influence  of  airfoil  leading  edge  radius  (as 
suggested  by  McCroskey^)  as  well  as  the  effect  of  different  static 
airfoil  pressure  gradients,  if  any,  for  the  case  of  large  amplitude 
motions . 


II.  DESCRIPTION  OF  THE  EXPERIMENT 

II. 1  MODEL  OSCILLATOR  SYSTEM 

An  investigation  of  the  effect  of  arbitrary  motions  on  the  dynamic 
stall  process  is  highly  dependent  on  the  development  of  a  suitable  model 
oscillation  mechanism  capable  of  offering  a  wide  range  of  motion  alterna¬ 
tives.  To  that  end  an  apparatus  capable  of  providing  simultaneous  motions 
of  an  airfoil  or  wing  model  in  both  pitch  and  translation  was  designed 
and  fabricated.  In  combination  with  the  measurement  techniques  described 
below,  the  concept  affords  a  high  degree  of  cycle-to-cycle  repeatability 
and  a  potential  to  generate  any  of  a  virtually  infinite  variety  of  simul¬ 
taneous  complex  motions  in  both  rotation  and  translation.  The  system 
was  specifically  designed  and  constructed  to  interface  with  the  U.S.  Air 
Force  Academy's  2  foot  (61  cm)  x  3  foot  (91  cm)  subsonic  wind  tunnel  - 
the  facility  in  which  these  studies  were  conducted.  A  detailed  descrip¬ 
tion  of  this  flow  facility  and  its  characteristics  are  provided  in 
reference  21. 

The  model  oscillation  system  is  functionally  composed  of  four  elements 
(1)  the  mechanical  assembly  or  "motion  transfer  apparatus,"  (2)  dual 
hybrid  servo  drive  units,  (3)  microcomputer  control  system,  and  (4)  the 
model  and  model  interface  collar.  These  elements  are  depicted  schemati¬ 
cally  in  Figure  II-l.  The  function  and  operation  of  each  of  these  elements 
will  be  described  briefly.  For  a  more  complete  discussion  of  the  oscilla¬ 
tion  system  and  its  characteristics,  the  reader  is  referred  to  reference 


II-2 


II. 1.1  Motion  Transfer  Apparatus 

The  motion  transfer  apparatus  is  a  mechanical  assembly  consisting 
of  various  drive  train  components  mounted  on  an  aluminum  frame  which 
transfers  mechanical  energy  from  the  drive  motors  to  the  wind  tunnel 
model.  During  operation,  the  unit  is  located  outside  of  and  below  the 
wind  tunnel  test  section  so  that  only  the  model  protrudes  into  the  air 
stream.  A  perspective  view  of  the  complete  apparatus  is  provided  in 
Figure  II-2,  as  viewed  from  the  downstream  quarter. 

The  drive  system  is  configured  in  a  dual  yoke  arrangement  which 
provides  uncoupled  motions  in  rotation  (pitching  motion,  in  this  case) 
and  a  single  degree  of  translation.  The  pitch  drive  assembly,  including 
the  pitch  motor,  is  "piggy-backed"  on  the  translational  drive  unit. 

Pitch  motion  is  generated  through  the  pitch  yoke  subassembly  depicted 
in  a  top  view  in  Figure  II-3.  The  mechanical  transfer  function  which 
relates  the  instantaneous  model  pitch  angle  (<* )  to  the  motor  shaft 
displacement  angle  (9j )  can  be  written: 


a 

max 


.  +2  li 

°  rl 


(2) 


The  entire  pitch  oscillation  mechanism,  including  model  and  drive 
motor,  is  mounted  directly  to  the  translational  drive  yoke  assembly. 

This  unit  is  shown  in  Figure  II-4  (top  view),  which  also  depicts  the 
relative  location  of  the  pitch  oscillator  components.  The  mechanical 
operation  of  the  translational  drive  train  is  functionally  similar  to 
that  of  the  pitch  assembly.  However,  it  was  not  employed  in  the  pitch- 
only  experiments  described  in  this  report,  and  will  not  be  further 
discussed.  Precision,  linear,  servo-mount  potentiometers  interfaced 
with  the  mechanical  drive  assemblies  serve  as  reference  position  trans¬ 
ducers.  Their  outputs  provide  model  angle  and  position  measurements 
with  less  than  0.052  offset  error  for  unidirectional  motions. 

II. 1.2  Hybrid  Servo  Drive  Units 

The  selection  of  a  suitable  drive  concept  to  power  the  two  oscillator 
stages  was  severely  impacted  by  two  demanding  performance  requirements, 
namely,  (1)  maximization  of  system  dynamic  response,  and  (2)  the 
selectable  variation  of  the  time  histories  of  the  motion  parameters, 
a(t)  and  x(t).  The  first  of  these  is  readily  attainable  through  the 
application  of  closed  loop  control  techniques.  The  latter  requirement, 
it  was  believed,  could  best  be  satisfied  through  the  use  of  real  time 
computer  control.  Ultimately,  both  requirements  were  met  by  employing 
both  of  these  ideas  and  utilizing  a  device  known  as  a  "servo-translator." 

The  operation  of  the  resultant  servo  system  utilizes  both  velocity 
and  position  feedback  to  achieve  effective  control.  Such  a  mechanism  is 
shown  schematically  in  Figure  II-5.  The  input  voltage  to  the  servo 
controller  is  a  "velocity  command"  which  is  instantaneously  and  continu¬ 
ously  compared  to  the  velocity  feedback  signal  provided  from  a  tachometer 


attached  to  the  drive  motor.  The  difference,  the  velocity-error  voltage, 
is  greatly  amplified  by  a  preamplifier  and  applied  to  a  power  amplifier 
to  command  motor  current.  The  current  generates  torque  which,  in  turn, 
alters  the  motor  speed  to  reduce  the  velocity  error  voltage.  Thus,  the 
servo  controller  is  always  trying  to  reduce  the  velocity  error  to  zero. 

A  separate  position  feedback  loop  controls  position  in  the  same 
manner  that  the  servo  controller  controls  speed.  A  position  command 
generated  by  an  external  source  is  compared  with  a  position  feedback 
signal  provided  by  a  linear  feedback  transducer  sensitive  to  position. 

Closed-loop  position  control  is  provided  by  an  input  stage  device 
termed  an  Absolute  Digital  Translator  Subtractor  (ADTS)*  which  accepts 
command  information  in  the  form  of  serial  data  pulses.  Pulse  information 
is  compared  with  the  feedback  signal  from  a  non-pulse  counting  type  of 
transducer  (termed  a  resolver)  to  close  the  loop.  This  design  increases 
noise  immunity  since  the  requirement  for  encoders  or  other  pulse-counting 
feedback  circuits  is  eliminated.  Errors  due  to  "accumulation"  associated 
with  these  techiques  are  also  eliminated. 

The  resultant  operating  system  is  controlled  in  real  time  by  the 
incoming  command  pulse  train.  Each  pulse  received  will  cause  the  motor 
to  move  one  increment  in  the  selected  direction,  and  the  rate  at  which 
the  pulses  are  received  determines  the  instantaneous  motor  speed.  The 
resolution  (magnitude  of  an  angular  increment)  is  regulated  in  the  design 
of  the  ADTS  stage  in  the  servo  translator.  Input  pulses  can  be  generated 
through  a  number  of  techniques  including  the  output  of  a  digital  computer 

*Patent  applied  for,  Control  Systems  Research  Corporation 
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II. 1.3  Microcomputer  Control  System 

The  choice  of  a  suitable  motion  control  computer  was  dictated  primarily 
by  the  nature  of  the  pulse  train  which  must  be  produced  as  an  input  for 
the  servo  translator.  To  generate  the  dual  pulse  trains  required  for 
simultaneous  operation  of  the  drive  systems,  a  Digital  Equipment  Corporation 
PDP  11/03  microcomputer  was  chosen  for  dedicated  motion  control.  The 
unit  was  configured  with  32  k-words  random  access  memory  (RAM)  which 
provided  operating  system  overhead  and  program  support. 

Control  of  the  two  motor  drive  systems  was  implemented  through  a 
single  digital  output  port  using  select  pins  for  the  various  pulse  and 
voltage  level  command  functions.  Pulses  were  generated  by  momentary 
"toggling"  of  the  voltage  at  the  appropriate  output  voltage  register  pin 
from  the  low  to  high  state.  Timing  of  the  pulses  for  each  control  channel 
was  accomplished  in  software  by  indexing  a  "pointer"  through  a  binary 
array.  In  this  process,  each  bit  is  examined  as  the  pointer  reaches  it, 
resulting  in  a  toggle  command  to  the  appropriate  output  port  if  a  "1"  is 
detected.  Bits  set  to  "0"  are  ignored.  Succeeding  bits  from  the  two 
arrays  which  determine  pitch  and  translational  motion  are  examined 
alternately.  The  pointer  indexing  rate  is  regulated  by  the  internal 
computer  instruction  clock,  and,  in  this  manner,  a  pulse  train  having  a 
prescribed  time  history  can  be  generated  programmatically.  Resulting 
motion  is  then  a  single  cycle  pattern  which  can  be  repeated  as  many 
times  as  required  by  relocating  the  pointer  to  the  top  of  the  arrays  and 
restarting  the  sequencing  program. 


II. 2  AIRFOIL  MODELS  AND  MODEL  INTERFACE 


Two  airfoil  shapes  were  employed  in  the  present  studies.  Initially, 
a  model  having  a  NACA  0012  cross-section  was  fabricated  with  rectangular 
planform  dimensions  of  14  inch  span  (35.6  cm)  and  6  inch  chord  (15.2 
cm).  The  model  was  constructed  from  three  solid  aluminum  pieces  and 
excavated  to  provide  interior  mounting  locations  for  nineteen  (19) 
miniature  piezo-resistive  pressure  transducers.  A  side  view  of  the 
airfoil  showing  the  pressure  measurement  locations  is  provided  in  Figure 
II-6. 

The  small  physical  size  of  the  model  precluded  flush-mounted  trans¬ 
ducers  and  necessitated  a  compromise  installation  method.  Short  pieces 
of  flexible  vinyl  tubing  were  used  to  interconnect  the  transducer  heads 
to  short  brass  tubes  which  extended  to  small  holes  in  the  model's  upper 
and  lower  surfaces.  Although  this  configuration  resulted  in  a  somewhat 
reduced  frequency  response,  the  upper  rolloff  limit  of  each  measurement 
channel  was  maintained  at  an  acceptable  level  ( fmax  1  1  kHz). 

To  examine  the  influence  of  airfoil  shape  on  deep  dynamic  stall 
behavior,  a  NACA  64]A012(13)  profile  was  also  selected  for  use  in  subsequent 
experiments.  This  profile  possesses  a  somewhat  smaller  leading  edge 
radius  than  the  NACA  0012.  McCroskey,  et  al?  have  suggested  this  as  a 
potentially  important  parameter  which  can  influence  the  evolution  of  the 
dynamic  stall  flowfield.  In  addition,  this  airfoil  possesses  a  longer 
favorable  pressure  gradient  region  over  the  forward  portion  of  the  surface 
characteristic  of  modern  aircraft  wing  sections.  For  example,  the  F-15 
aircraft  employs  the  12%  thickness-to-chord  ratio  version  of  this  airfoil 
shape  in  the  wing  root  region. 
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Construction  of  the  second  airfoil  was  similar  to  that  of  the  original 
except  for  the  physical  size  -  14  inch  span  (35.6  cm)  with  an  8  inch 
chord  (20.3  cm).  The  slightly  larger  chord  of  the  second  model  was 
chosen  to  increase  the  effective  spatial  resolution  and  improve  the  ease 
of  transducer  installation.  It  was  designed  to  accommodate  20  internally 
mounted  pressure  transducers  as  shown  in  Figure  11^7. 

Note  that  neither  model  was  constructed  to  fully  span  the  wind 
tunnel  test  section.  Limits  on  the  practical  spanwise  length  were 
imposed  by  inertial  considerations  which  influence  realizable  model 
accelerations  during  unsteady  motions.  To  approximate  true  airfoil 
characteristics  and  maximize  flow  two-dimensionality ,  plexiglass  endplates 
were  constructed  for  both  models.  These  nearly  rectangular,  spanwise 
adjustable  pieces  were  configured  to  project  at  least  0.20  chord  lengths 
upstream  and  downstream  of  the  leading  and  trailing  edges,  respectively. 
They  also  extended  approximately  one-half  of  a  chordlength  in  the  cross¬ 
stream  direction  from  the  upper  and  lower  airfoil  surfaces.  When  posi¬ 
tioned  for  the  tests  described  in  subsequent  sections,  they  were  placed 
equidistant  from  the  pressure  tap  centerline  to  give  an  effective  geometric 
aspect  ratio  of  1.3  for  both  airfoil  models.  Care  was  exercised  to 
place  the  lower  endplate  (nearest  the  wind  tunnel  floor)  out  of  the  test 
section  wall  boundary  layer. 

Both  models  were  designed  to  accommodate  an  interface  collar  which 
could  be  coupled  directly  to  the  pitch  drive  assembly.  The  pitch  drive 
spur  gear  mentioned  in  Section  II. 1  is  permanently  attached  to  the  hollow 
model  support  shaft  which  is,  in  turn,  fastened  to  the  lower  portion  of 
the  aluminum  mainframe  by  two  alignment  bearings.  The  model  is  coupled 
to  the  shaft  through  a  'tang'  which  protrudes  from  the  base  of  the  model. 


The  tang  was  fabricated  to  be  inserted  into  the  upper  opening  of  the 
support  shaft  and  held  in  place  with  recessed  set  screws.  The  tang 
posesses  an  oblong  cross-section  and  is  hollow  to  accomodate  transducer 
wires  passing  from  the  model's  interior  to  outside  electronic  signal 
conditioning  instrumentation.  A  sketch  of  the  model  mount  arrangement 
is  provided  in  Figure  II-8. 

Due  to  model  size  and  shape  limitations,  the  pitch  rotation  points 
of  the  two  airfoils  are  slightly  different.  The  NACA  0012  model  rotates 
about  the  31.7  per  cent  chord  point  while  the  NACA  64iA012(13)  model 
rotates  at  the  37.5  per  cent  chord  location.  Previous  studies  have 
indicated  that  the  location  of  the  pitching  axis  has  little  effect  on 
the  peak  loading  value  but  can  influence  the  angle  at  which  stall  actually 
occurs. ^  Comparisons,  in  this  case,  should  be  made  with  caution. 

II. 3  PRESSURE  MEASUREMENT  SYSTEM 

Unsteady  pressures  were  measured  by  miniature,  dynamic  piezo-resistive 
transducers  possessing  a  high  frequency  response.  For  the  current  experi¬ 
ments,  Endevco  Model  8507-2  transducers  specified  over  nearly  a  100  kHz 
bandwidth  and  nominally  yielding  120  mv/psi  sensitivity  at  10  volts 
excitation  were  employed.  The  specified  sensitivity  of  these  trans¬ 
ducers  to  acceleration  is  minimal  at  0.0005  psi/g,  maximum.  The  devices 
are  also  temperature  compensated  to  provide  minimal  drift  from  this 
sensitive  parameter  during  operational  use.  The  short  piece  of  flexible 
tubing  connecting  each  transducer  with  the  appropriate  surface  measure¬ 
ment  port  reduced  the  overall  upper  frequency  response  limit  to  levels 
no  lower  than  1  kHz  in  the  worst  case. 


To  provide  a  stable  reference  for  the  unsteady  flow  experiments, 
the  "reference"  side  of  each  of  the  transducers  was  connected  to  atmos¬ 
pheric  pressure  in  the  room  exterior  to  the  wind  tunnel.  This  arrange¬ 
ment  was  judged  preferable  to  alternative  schemes  using  in-tunnel 
reference  sources  which  might  otherwise  introduce  errors  due  to  unsteady 
phase  lag  variations.  Steady  flow  total  and  static  pressures  in  the  test 
section  were  measured  using  Statham  Model  PM5TC  transducer  units.  An 
interconnection  diagram  of  the  pressure  measurement  system  is  provided 
in  Figure  11-9.  Using  this  system,  the  transducer  outputs  sampled  were 
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The  unsteady  pressure  coefficients  were  computed  using  the  following 
relationship: 
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where  the  "o"  subscript  denotes  reference  steady  flow  conditions. 

Each  of  the  dynamic  pressure  transducers  was  excited  by  a  separate, 
regulator-isolated,  DC  power  supply  at  10  volts.  A  transducer  bridge 
balance  circuit  provided  coarse  and  fine  adjustments  of  the  desired 
reference  level  of  the  individual  output  voltages.  The  low-level  output 
signals  were  directed  to  preamplifier  units  located  at  the  extremities 


of  the  transducer  lead  wires  to  minimize  noise  interference  effects. 
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The  preamplified  signals  were  then  channelled  to  Dynamics  Model  7521B 
variable  gain  amplifier  units  to  optimize  the  signal  level  to  a  value 
compatible  with  the  data  acquisition  system.  A  complete  schematic  of 
the  electronic  signal  conditioning  system  is  provided  in  Figure  11-10. 

II. 4  DATA  ACQUISITION  SYSTEM 

Although  the  microcomputer  control  system  described  earlier  was 
adequate  to  achieve  motion  control,  its  limited  memory  and  lack  of 
storage  peripherals  severely  restricted  its  use  in  program  development. 
Therefore,  a  larger  minicomputer  (DEC  Model  PDP  11/45)  was  employed  for 
overall  experimental  control  and  data  acquisition.  The  central  processor 
was  configured  with  128k  words  of  memory  (16  bits  per  word)  and  20  mega¬ 
words  of  disk  storage  capacity.  [Of  this  capacity,  a  maximum  of  32k-words 
of  memory  were  available  to  a  single  user  under  the  RSX  11-M  operating 
system. ]  Additional  mass  storage  was  accessible  through  a  network  link 
(DECNET)  to  other  PDP  11  computing  systems,  as  well  as  a  magnetic  tape 
drive  and  other  peripherals. 

This  system  was  equipped  with  a  laboratory  peripheral  system  (DEC 
Model  LPS-11)  having  a  32  analog  channel  input  capacity,  as  well  as 
three  16-bit  digital  input-output  ports.  The  LPS-11  module  was  equipped 
with  a  high-speed  real  time  clock,  multiple  Schmitt  triggers,  relays  and 
other  hardware  interface  subsystems  designed  for  use  in  laboratory 
applications. 

A  capability  was  developed  to  interface  the  larger  machine  with  the 
smaller  one  for  transfering  task  images  using  a  "down-line  load"  technique. 
The  simple  network  allowed  for  program  development  and  editing  in  the 
larger  minicomputer  system.  A  more  complete  discussion  of  the  elements 
of  motion  control  can  be  found  in  reference  22. 
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II. 5  MEASUREMENT  PROCEDURES 

Unsteady  pressure  measurements  during  actual  model  tests  were 
preceded  by  an  involved  calibration  and  alignment  sequence  designed  to 
minimize  sources  of  potential  measurement  error.  This  procedure  began 
with  the  careful  adjustment  and  alignment  of  all  mechanical  components 
to  minimize  friction  and  wear  and  to  maximize  cycle-to-cycle  repeat¬ 
ability  of  the  model  motion  waveform.  Several  reference  variables 
including  the  initial  angle  of  attack,  reference  steady  freestream 
velocity  and  other  reference  pressures  required  in  the  experiment  were 
initialized  and  measured  prior  to  each  test  sequence.  The  maximum 
incidence  angle  was  also  independently  verified.  Immediately  prior  to 
the  "flow  on"  condition,  adjustments  of  the  zero  set  point  (reference 
condition)  for  each  of  the  dynamic  transducers  were  accomplished.  These 
reference  values  were  rechecked  immediately  after  concluding  the  test 
to  determine  the  output  drift  of  all  data  channels.  A  variation  of  less 
than  1%  of  the  maximum  measured  level  for  any  channel  was  required  to 
validate  the  data. 

During  the  experiments,  the  typical  response  amplitudes  of  the 
dynamic  pressure  transducers  reached  a  small  fraction  of  their  rated 
value  (+  2  psid)  rendering  the  manufacturer's  calibration  constants 
suspect.  As  a  result,  individual  "system"  calibrations  were  conducted 
for  each  of  these  units  on  a  regular  basis  to  minimize  calibration 
errors  as  well  as  to  identify  defective  instrumentation  which  may  have 
been  damaged  during  prior  experiments.  With  one  exception  (failure  of 
transducer  #19  -  lower  surface,  near  trailing  edge  -  in  the  NACA  0012 
airfoil  during  an  early  test  series),  no  calibration  drift  variations 
greater  than  3%  were  detected  during  the  experiments. 
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In  addition  to  the  involved  pressure  and  velocity  measurement 
calibrations,  extreme  c are  was  also  exercised  in  the  assessment  of  model 
motion  parameters.  Motion  repeatability  was  assessed  both  prior  to  and 
following  each  of  the  tests.  In  all  cases  where  appropriate  care  in 
alignment  had  been  exercised,  motion  variations  were  not  discernable  on 
an  oscillograph  trace  (less  than  0.5%  variation  between  the  initial  and 
final  curves.) 

The  actual  acquisition  of  data  involved  repetitive  sampling  of  each 
of  the  pressure  channels  over  25  cycles  to  obtain  an  ensemble  average 
before  proceeding  to  the  next  adjacent  transducer.  This  process  was 
repeated  around  the  airfoil  until  all  transducers  had  provided  the 
requisite  output.  The  model  motion  time  history  was  also  stored  as  part 
of  the  data  file. 

II. 6  DATA  REDUCTION  AND  ANALYSIS 

Data  were  analyzed  using  several  programs  developed  for  the  PDP 
11/45  data  acquisition  system.  Ensemble-averaged,  unsteady  pressure 
time  histories  for  each  transducer  were  stored  as  elements  of  a  file 
representing  a  given,  unique  test  condition.  These  results  were  subse¬ 
quently  converted  to  time-frozen,  pressure  distributions  representative 
of  fixed  temporal  phase  points  during  the  motion  cycle. 

The  integration  of  the  surface  pressure  data  was  accomplished  with 
an  algorithm  based  on  the  trapezoidal  rule  and  appropriate  coordinate 
transformations  to  convert  normal  and  chord  forces  to  lift  and  drag 
forces.  This  process  incorporated  local  values  of  the  model  surface 
slope,  as  well  as  the  instantaneous  chord  line  orientation  to  determine 
the  appropriate  transformation  angles.  The  integration  of  these  lift 
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coefficient  data  to  calculate  the  life-impulse  function  (see  Section  V) 
was  also  accomplished  using  similar  numerical  procedures. 

No  attempt  was  made  to  apply  wind  tunnel  corrections  to  any  of  the 
unsteady  pressure  data.  The  original  choice  of  relatively  small  airfoil 
models  was  prompted  by  the  desire  to  minimize  errors  due  to  wall  effects. 

II. 7  SOURCES  OF  EXPERIMENTAL  ERROR 

The  influence  of  the  wind  tunnel  wall  proximity  on  the  measurements 
is  undoubtedly  a  factor  of  concern  for  any  assessment  of  dynamic  forces, 
especially  when  the  motion  amplitudes  are  large  as  in  the  present  case. 
Tunnel  blockage,  lift  effect  and  other  conventional  steady  flow  corrections 
are  expected  to  be  augmented  by  unsteady  counterparts  whose  magnitudes  are 
unknown  and  largely  unpredictable  given  the  current  state  of  unsteady  flow 
measurement  technology.  Correction  factors  would  be  expected  to  be  func¬ 
tions  of  time  (or  motion  phase  angle). 

The  periodic  (albeit,  anharmonic)  character  of  these  motions  introduces 
another  possible  source  of  "contamination,"  namely,  the  influence  of 
flow  phenomena  which  are  either  multiples  or  subharmonics  of  the  oscilla¬ 
tion  period.  Related  acoustic  effects  which  might  affect  the  behavior 
of  dynamic  Lagrangian  vortices  could  conceivably  introduce  sizeable 
errors.  It  is  believed  that  these  error  sources  were  minimized  by  verifying 
that  the  convective  tunnel  time  constant  was  sufficiently  large  compared 
to  the  period  of  model  motion  and  also  insuring  that  it  was  not  an  integral 
multiple  of  that  value.  No  perceptible  acoustic  interference  was  detected 
during  the  conduct  of  any  of  the  measurements.  It  is  believed  that  the 
very  energetic  character  of  the  separation  vortices  will  overwhelm  any 
weak  excitation  provided  by  experimentally  related  acoustic  field  effects 


at  least  during  the  important,  high  lift  portions  of  the  cycle  although 
this  assertion  is  speculative. 

The  selection  of  the  record  size  necessary  for  an  accurate  representa 
tion  of  the  phase-locked  average  flow  behavior  was  another  issue  of 
concern.  A  related  issue  was  the  accompanying  cycle  deviation  from  this 
average.  A  previous  investigation^  has  suggested  that  a  50  record 
ensemble  is  suitable  in  tests  of  this  type.  Due  to  the  need  to  minimize 
electronic  drift  and  given  the  lengthy  measurement  times  anticipated 
(several  hours),  it  was  decided  to  compromise  at  a  total  of  25  records 
for  each  transducer  to  calculate  the  phase  averaged  pressures.  The 
choice  of  this  value  was  the  result  of  a  series  of  preliminary  tests 
under  several  model  oscillation  conditions  which  examined  1,  5,  10,  20, 

25,  50  and  100  record  averages.  Although  the  25  record  average  is  a 
compromise,  worst  case  errors  resulting  from  this  limitation  were  esti¬ 
mated  not  to  exceed  10%  even  under  turbulent  flow  conditions. 

A  final  factor  which  bears  on  the  validity  and  extendability  of 
these  results  is  that  of  the  Reynolds  number  regime  in  which  the  experi¬ 
ments  were  conducted.  Chord  Reynolds  numbers  in  the  range  of  the  present 
measurements  and  associated  boundary  layer  development  are  known  to  lie 
within  what  is  considered  the  "transitional"  regime  for  typical  steady 
flow  conditions.  However,  it  is  believed  that  the  violent  large  amplitude 
motions  employed  in  the  deep  dynamic  stall  experiments  described  here 
will  readily  dominate  local  shear  layer  instabilities  immediately  prior 
to  stall  and  afterward  when  energetic  separation  vortices  are  present. 
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III.  AIRFOIL  PRESSURE  DISTRIBUTIONS 


Dynamic  pressure  measurements  have  been  obtained  for  both  airfoil 
models  over  a  wide  parametric  range  of  model  motion  and  reference  flow 
conditions.  Displayed  in  the  form  of  instantaneous  distributions  over 
the  airfoil  surface,  these  data  are  useful  not  only  for  quantifying 
airfoil  performance  as  a  function  of  time,  but  also  assist  in  the  detec¬ 
tion  of  anamolous  flow  behavior  characteristic  of  dynamic  stall  phenomena. 
Phase  locked,  ensemble -averaged  values  of  the  pressure  coefficient  are 
presented  as  functions  of  the  chordwise  position  coordinate  in  the  displays 
discussed  below. 

III.l  REFERENCE  CONDITIONS  -  STEADY  FLOW 

Before  conducting  surface  pressure  measurements  in  the  unsteady 
flow  environment,  data  were  obtained  Cor  both  airfoils  at  prescribed, 
fixed  incidence  angles  under  steady  freestream  conditions.  Some  of 
these  results  are  displayed  in  Figures  III-l  and  III-2  for  representative 
freestream  conditions  and  for  several  values  of  the  incidence  angle. 

These  data  were  obtained  to  further  confirm  the  accuracy  and  test  the 
stability  of  the  dynamic  pressure  transducer  measurement  system  and  to 
serve  as  a  reference  data  source  for  comparison  with  unsteady  flow 
results . 

In  these  figures  and  in  other  computer-generated  graphics  displayed 
in  this  report,  experimental  data  points  are  not  explicitly  identified 
but  have  been  connected  with  straight  line  segments  in  the  graphics 
software  program.  Actual  data  points  can  be  found  by  observing  those 
locations  where  the  interconnecting  straight  line  segments  abruptly 
change  slope  or,  alternatively,  noting  the  airfoil  pressure  measurement 
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locations.  Also  note  that  the  polarity  of  the  ordinate  is  reversed 
between  each  of  the  two  figures.  All  remaining  pressure  distributions 
discussed  in  this  report  are  displayed  with  positive  pressure  increasing 
upward  on  the  page  (in  the  sense  of  Figure  III-2). 

Although  results  for  only  a  single  freestream  velocity  are  depicted 
for  each  airfoil,  data  over  a  range  of  freestream  speeds  were  observed 
to  exhibit  identical  trends  with  changing  incidence  angle.  These  distri¬ 
butions  were  found  to  be  remarkably  insensitive  to  the  chord  Reynolds 
number  over  the  range  examined,  except  for  results  obtained  near  the 
static  stall  angle. 

The  large  rearward  extent  of  the  adverse  pressure  gradient  on  the 
upper  airfoil  surface  at  moderate  to  high  incidence  angles  is  not 
unexpected  for  the  classes  of  airfoils  employed.  The  presence  of  a 
separation  "bubble"  over  the  first  quarter  chord  is  also  anticipated 
at  higher  incidence  angles.  The  region  of  favorable  pressure  gradient 
on  the  upper  (suction)  surface  is  more  extensive  for  the  64-series 
airfoil  at  lower  angles  of  attack. 

Results  of  the  NACA  0012  measurements  of  Figure  III-l  have  been 
compared  with  Navier-Stokes  computations  by  Hegna^S,  who  obtained 
numerical  results  for  identical  reference  flow  conditions.  Good  quali¬ 
tative  and  quantitative  agreement  was  achieved  with  a  slight  discrepancy 
in  the  trend  of  the  distributions  toward  the  airfoil  trailing  edge.  The 
overall  magnitude  of  disagreement  has  been  approximated  to  be  within  the 
range  of  5  to  8  per  cent  for  all  cases  where  a  comparison  is  available. 

Integral  results  associated  with  these  distributions  are  discussed 


in  Section  IV. 


III. 2  UNSTEADY  TEST  CONDITIONS 


A  complete  matrix  of  unsteady  test  conditions  employed  in  this 
series  of  experiments  is  provided  for  both  airfoils  in  Tables  III-l  and 
III-2  in  the  Appendix.  The  model  motion  parameters  listed  are  based  on 
measured  performance  as  determined  from  data  from  pitch  and  translation 
position  transducers.  The  test  matrix  was  designed  to  maximize  the 
range  of  the  test  parameters  which  could  be  varied  in  assessing  airfoil 
performance  throughout  the  separation  regime.  Although  there  is  some 
overlap  in  the  motion  range  for  the  two  airfoils  (nearly  identical  test 
conditions),  there  are  also  a  significant  number  of  cases  for  which  no 
comparison  is  available. 

A  graphical  depiction  of  various  motion  waveform  shapes  is  presented 
in  Figure  III-3  to  assist  in  the  interpretation  of  reference  data  provided 
in  the  Tables.  The  use  of  anharmonic  motions  greatly  facilitated  the 
attainment  of  large  amplitudes  at  comparatively  high  instantaneous  pitch 
rates,  a  major  objective  of  this  investigation.  Although  a  wide  variety 
of  motion  histories  have  been  examined,  this  report  addresses  only  those 
cases  involving  pure  pitching  motions  at  constant  rate,  i.e. ,  «  =  constant. 
The  remaining  motions,  including  combined  two-degree-of- freedom  experiments, 
will  be  addressed  in  a  subsequent  technical  report. 

Large  amplitude,  constant  rate,  pitching  motions  such  as  those 
characterized  by  case  A  in  Figure  III-3  were  the  subject  of  study  for 
several  reasons.  These  motions  are  capable  of  providing  a  direct  correla¬ 
tion  between  what  is  believed  to  be  the  most  significant  motion  variable, 
pitch  rate  (°0,  and  leading  edge  separation  behavior.  Unlike  the  many 
previous  dynamic  stall  studies  involving  harmonic  pitch  oscillations, 
the  development  of  the  leading  edge  vortex  is  not  'forced'  by  a  decrease 
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in  airfoil  pitch  rate  (“  "*■  0)  near  a  relatively  low  maximum  incidence 
angle.  The  stall  process  develops  as  the  motion  continues  unaltered 
well  beyond  the  onset  of  separation.  The  occurrence  of  leading  edge 
separation  and  stall  vortex  development  prior  to  reaching  the  maximum 
incidence  then  represents  a  unique  physical  situation  of  special 
interest  to  this  investigation. 

In  contrast,  the  idea  of  "trapping"  a  separation  vortex  to  achieve 
some  measure  of  lift  enhancement  requires  that  the  airfoil  "respond"  to 
the  vortex  after  it  has  formed.  For  this  purpose,  it  was  thought  that 
the  sudden  discontinuation  of  the  nose  up  pitching  motion  at  some 
optimum  intermediate  angle  might  delay  the  passage  of  the  vortex  in  the 
blockage  region  behind  the  high  incidence  configuration. 

The  pursuit  of  lift  enhancement  under  post-stall  conditions  was 
also  the  motivation  for  studying  motions  such  as  those  shown  in  case  B 
(constant  rate  pitch  motion  to  maximum  incidence  angle  followed  by  an 
abrupt  decrease  in  incidence).  In  this  case,  the  sudden  reduction  in 
pitch  angle  was  accomplished  to  strengthen  (or  energize)  the  shed  stall 
vortex  by  increasing  the  local  shear  stress  near  the  airfoil  leading 
edge. 

In  all  cases,  the  cycle  duration  was  prolonged  to  allow  for  a  relaxa¬ 
tion  of  induced  unsteady  effects  following  the  completion  of  motion. 

For  example,  subsequent  to  a  constant  pitch  rate  motion  (type  A),  the 
airfoil  was  kept  in  the  maximum  incidence  position  until  static  stall 
conditions  were  achieved.  This  constraint  insured  that  cycle-to-cycle 
coupling  effects  induced  by  neighboring  vortex  interactions  were 


eliminated. 


To  insure  a  complete  relaxation  to  steady  flow,  stall  conditions 
during  each  cycle  of  motion,  the  delay  time  for  airfoil  position  reset; 
that  is,  the  time  between  the  initiation  of  motion  and  the  reset  command 
was  calculated  from  the  following  relationship: 

t  >!?£*+  3c  (7) 

d  —  a  U 

OO 

A  range  of  pitch  rates  over  nearly  three  orders  of  magnitude  was 
employed  to  study  constant  rate  motion  for  the  NACA  0012  airfoil. 

Maximum  incidence  angles  ranging  from  nearly  20  degrees  to  over  60  degrees 
were  considered.  The  low  end  of  the  rate  and  amplitude  ranges  extended 
from  conditions  studied  by  previous  dynamic  stall  investigators  to  an 
opposite  extreme  in  which  incidence  angles  well  in  excess  of  the  induced 
separation  orientation  could  be  achieved  for  even  the  highest  pitch 
rates  generated.  A  generally  higher  range  of  pitch  rates  was  employed 
for  the  NACA  64iA012(13)  airfoil  experiments. 

III. 3  SHAPE  OF  THE  DISTRIBUTIONS 

Despite  the  variations  of  motion  type  and  the  range  of  motion 
parameters  employed  in  the  experiments,  many  of  the  qualitative  features 
associated  with  the  shapes  of  the  surface  pressure  distributions  as  they 
evolve  in  time  are  similar  over  the  entire  measurement  range.  To  examine 
this  behavior,  consider  the  typical  historical  record  depicted  in  Figure 
III-4.  This  example  illustrates  a  case  where  the  apparent  separation 
angle  is  exceeded  by  the  maximum  amplitude  of  model  motion.  At  values 
less  than  the  static  stall  angle  (a),  the  shape  of  the  pressure  distribu¬ 
tion  is  similar  to  that  expected  for  steady  conditions  except  for  an 
apparent  time  lag  in  the  pressure  coefficient  magnitudes.  As  the  inci- 
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dence  angle  increases  to  and  beyond  the  static  stall  value,  an  "attached 
flow"  distribution  persists  with  suction  pressure  magnitudes  in  excess 
of  equivalent  steady  flow  values. 

As  the  incidence  angle  continues  to  increase  (b),  the  suction  peak 
is  observed  to  "break  away"  from  the  leading  edge  region  on  the  airfoil 
upper  surface  and  move  downstream  along  the  surface.  This  behavior  is 
highly  suggestive  of  the  dynamic  stall  vortex  observed  in  previous 
airfoil  experiments  since  the  evolution  of  the  distribution  is  identical 
and  in  the  same  relative  range  of  incidence  angles  for  comparable  pitch 
rates.  The  increased  amplitude  of  the  suction  peak  also  suggests  that 
the  "vortex"  is  energized  as  the  incidence  angle  continues  to  increase 
over  a  range  determined  by  the  model  rotation  rate  (<*). 

At  incidence  angles  well  in  excess  of  the  static  stall  value  (c), 
the  pressure  peak  begins  to  decrease  in  magnitude  and  broaden  over  the 
surface  as  it  continues  to  move  rearward.  A  possible  interpretation  of 
this  behavior  is  that  the  dynamic  stall  vortex  begins  to  lose  energy  as 
it  continues  to  translate  into  the  blocked  "wake  region"  created  by  the 
large  pitch  angle  of  the  airfoil.  As  this  rotating  fluid  loses  energy, 
the  vortex  structure  becomes  more  diffuse  as  it  ultimately  evolves  into 
a  separated  free  shear  layer  emanating  from  the  airfoil  leading  edge. 

As  the  motion  continues  further  (or,  in  some  high  rate  cases,  after 
the  motion  is  halted),  the  flow  field  relaxes  to  a  state  which  is  essen¬ 
tially  that  of  a  bluff  body  wake  (d).  The  pressure  distributions  over 
the  upper  and  lower  surfaces  flatten  considerably  and  attain  levels  of 
nearly  plus/minus  unity,  respectively. 

Experimental  results  confirm  that  the  extent  of  these  behavioral 
trends  is  a  strong  function  of  the  model  motion  parameters,  especially 


the  pitch  rate  and  the  pitch  amplitude.  Increased  pitch  rate  tends  to 
delay  the  onset  of  leading  edge  unsteady  separation  and,  in  general, 
extend  the  apparent  attached  flow  behavior  to  higher  incidence  angles. 

As  suggested  earlier,  the  maximum  incidence  angle  can  have  a  major  impact 
on  the  shape  of  the  distributions  especially  if  it  does  not  exceed  the 
value  at  which  the  separation  would  otherwise  occur.  These  amplitude 
effects  are  addressed  more  completely  in  sections  IV  and  V. 

Representative  distributions  for  both  airfoils  are  displayed  in 
Figures  III-5  through  III-10.  These  data  facilitate  a  comparison  of 
separated  flow  performance  for  both  airfoils  at  three  comparable  levels 
of  reduced  pitch  rate.  In  all  cases,  the  angle  for  which  unsteady 
separation  occurs  is  well  exceeded  by  the  amplitude  of  motion,  amax>  of 
approximately  60  degrees.  The  shapes  of  the  distributions  for  both 
airfoil  models  are  found  to  exhibit  the  same  general  qualitative  features. 
However,  the  suction  peak  associated  with  the  NACA  64-series  airfoil  was 
observed  to  be  steeper  and  more  pronounced  at  high  pitch  rates  and  was 
also  observed  to  "detach"  from  the  leading  edge  at  high  incidence  angles. 

The  shape  of  the  NACA  0012  distribution  did  not  exhibit  the  same  accentuated 
behavior  under  similar  conditions. 

Additional  similarities  in  the  pressure  time  histories  become  apparent 
through  an  examination  of  the  amplitude  of  the  minimum  pressure  coeffi¬ 
cient  associated  with  the  suction  peak  as  it  evolves  with  increasing 
incidence  angle.  These  results  are  displayed  in  Figures  III-ll  through 
III-13.  Although  some  differences  in  amplitude  levels  are  apparent,  the 
results  are  observed  to  be  in  good  agreement  when  compared  to  the  varia¬ 
tion  in  dimensionless  pitch  rate. 
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IV.  INTEGRAL  MEASUREMENT  RESULTS 


The  pressure  distributions  described  in  the  preceding  section  have 
been  integrated  to  provide  load  and  moment  time  histories  for  each  of 
the  test  cases  examined.  As  indicated  in  Section  II,  numerical  integra¬ 
tion  programs  provided  values  of  airfoil  normal  force  coefficient  (Cn), 
lift  coefficient  (C^),  moment  coefficient  about  the  quarter  chord  (C„,), 
and  pressure  drag  coefficient  (C<j)  aa  functions  of  time  and,  implicitly, 
incidence  angle.  Available  computer  graphics  capabilities  provided 
cross  pic  -  of  these  variables,  with  time  and  incidence  angle  as  potential 
parameters . 

Graphs  of  lift  coefficient  versus  incidence  angle  for  steady  free- 
stream  and  static  modal  conditions  are  displayed  for  both  airfoils  in 
Figures  IV-1  and  IV-2.  Maximum  values  of  the  static  lift  coefficient 
are  alao  shown.  NACA  0012  data  are  plotted  with  results  of  a  Navier- 
Stokes  computation2^  and  are  observed  to  be  in  good  agreement.  The  low 
value  of  the  initial  lift  curve  slope  as  well  as  the  nonlinear  behavior 
at  higher  incidence  angles  are  believed  to  be  attributable  to  the  test 
Reynolds  number  range.  Note,  however,  that  the  variation  in  lift 
characteristics  over  the  range  of  Reynolds  numbers  studied  here  (factor 
of  2)  is  small  compared  to  the  departure  from  ideal  airfoil  behavior. 

Unsteady  lift,  moment,  and  drag  data  for  both  airfoils  are  displayed 
in  Figures  IV-3  through  IV-12  for  several  sets  of  reference  conditions. 
These  data  sets  were  selected  to  afford  a  basis  for  comparing  performance 
of  the  two  airfoils  at  nearly  identical  pitch  rates  and  amplitudes. 

Note  that  the  lift  coefficient  is  displayed  instead  of  the  more  tradi¬ 
tional  normal  force  coefficient,  since  the  primary  objective  of  this 


report  is  an  assessment  of  aerodynamic  lifting  performance.  Normal  and 
chord  force  coefficients  have  been  converted  to  lift  and  drag  coeffi¬ 
cients  through  a  cartesian  coordinate  transformation  which  accounts  for 
the  instantaneous  airfoil  orientation. 

Lift  stall  characteristics  are  observed  to  be  comparable  for  the 
two  airfoils  at  low  and  moderate  pitch  rates  (k  «*  0.045  and  0.09, 
respectively)  with  large  amplitude  displacements.  However,  the  NACA 
0012  airfoil  is  observed  to  provide  a  slightly  less  severe  moment  stall 
(large  nose  down  moment)  in  this  range.  When  the  maximum  amplitude  of 
motion  is  less  than  or  comparable  to  the  incidence  angle  at  which  unsteady 
leading  edge  separation  would  otherwise  occur  (Figures  IV-9  through 
1V-12),  the  moment  stall  associated  with  both  airfoils  is  observed  to  be 
comparable  and  severe  at  the  maximum  incidence  angle.  An  abrupt  decrease 
in  moment  coefficient  is  also  evident  in  the  plot  of  C,,,  with  time.  The 
overall  shape  of  the  lift  curves  is  similar  for  the  two  airfoils  under 
any  comparable  set  of  flow  and  motion  conditions. 

The  pressure  drag  force  is  also  important  in  assessing  overall 
airfoil  performance.  Drag  coefficients  of  relatively  large  magnitude 
are  observed  in  conjunction  with  dynamic  stall;  a  fact  which  is  physically 
consistent  with  a  notion  that  energetic  unsteady  separation  results  in  a 
tightly  wound  vortex-like  structure  during  the  initial  stages  of  the 
process.  If  the  airfoil  were  performing  in  a  free  flight  environment, 
its  large  increasing  drag  force  would  result  in  an  axial  deceleration  in 
the  freestream  direction.  Alternatively,  it  can  be  hypothesized  that 
the  transient  lift  associated  with  the  dynamic  stall  vortex  is  provided 
through  an  extraction  of  freestream  kinetic  energy  due  to  the  airfoil 
motion.  The  mechanics  of  this  process  may  be  more  fully  understood 
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by  considering  airfoil  pitching  motions  which  are  coupled  to  time 
dependent  behavior  of  the  freestream.  A  discussion  of  some  preliminary 
results  which  address  this  issue  will  be  the  subject  of  a  later  report. 

The  shape  of  the  drag  coefficient  curves  are  similar  for  the  two 
airfoils  with  one  notable  exception.  Initially,  the  drag  increases  less 
rapidly  for  the  NACA  64]A012(13)  airfoil  as  the  incidence  angle  rises. 
However,  this  is  followed  by  a  rapid  increase  in  the  drag  coefficient  to 
the  maximum  value.  The  rate  of  drag  increase  of  the  NACA  0012  is  more 
uniform  over  the  range  of  angles  of  attack. 

The  occurrence  of  separation  and  stall  prior  to  reaching  the  maximum 
incidence  of  motion  is  again  evident  from  the  time  history  data  shown  in 
Figure  IV-13.  These  results  suggest  that  the  unsteady  separation  "vortex" 
is  shed  at  some  lower  angle  which  coincides  with  the  formation  and  progres¬ 
sion  of  the  pressure  suction  peak.  Force-time  curves  of  this  type  are 
useful  in  assessing  overall  airfoil  performance  as  discussed  in  the  next 
chapter. 

An  interesting  feature  apparent  in  the  lift  coefficient  versus  time 
curves  is  that  of  a  "notch"-like  indentation  which  is  noticeably  larger 
than  the  more  random  smaller  amplitude  fluctuations.  The  notch  is 
depicted  schematically  in  Figure  IV-14.  It  occurs  as  a  depression  in 
the  lift  coefficient  between  two  local  maxima  and  may  well  be  caused  by 
the  appearance  of  a  secondary  vortex  structure  or  possibly  other  vortex 
behavior  which  must  be  documented  through  flowfield  measurements.  Using 
the  notation  illustrated  in  Figure  IV-14,  the  notch  location  (time 
difference  between  its  occurrence  and  that  of  the  maximum  lift)  is 
plotted  as  a  function  of  dimensionless  frequency  in  Figure  IV-15.  Notch 
location  is  observed  to  be  only  a  weak  function  of  frequency  but  is 
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different  for  the  two  airfoils.  For  cases  of  high  maximum  incidence 
(d^ax  *  60  degrees)  and  low  pitch  rate,  the  notch  is  observed  to  occur 
in  advance  of  the  time  that  the  airfoil  reaches  its  maximum  lift  coeffi¬ 
cient.  With  increasing  pitch  rate,  the  notch  is  observed  to  appear 
later  in  the  cycle  and  eventually  will  occur  subsequent  to  the  maximum 
lift  point.  This  result  is  especially  apparent  for  the  case  of  the  HACA 
64]A012(13)  airfoil.  The  rate  of  shift  in  notch  location  for  the  two 
airfoils  is  roughly  comparable  for  high  levels  of  maximum  incidence  in 
the  region  of  overlap.  The  notch  phenomenon  was  also  observed  for  cases 
of  lower  maximum  incidence  angle  (amax  “  35-40  degrees).  However,  in 
those  cases,  it  was  always  observed  to  occur  after  the  attainment  of 
maximum  lift.  The  repeatability  of  this  behavior  and  persistence  of  the 
notch-like  depression  suggests  that  the  responsible  viscous-inviscid 
interaction  occurs  regularly  on  airfoils  undergoing  deep  dynamic  stall. 
This  behavior  is  recommended  as  the  subject  of  future  study. 


STEADY  FLOW  LIFT  COEFFICIENTS,  NACA  0012  AIRFOIL 
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FIGURE  1V-4 

FORCE  AND  MOMENT  CHARACTERISTICS,  NACA  64,A012{13)  ,  k  =  0.045,amax  =  58°,  UL=  32.1  fps. 


FIGURE  IV -7 

FORCE  AND  MOMENT  CHARACTERISTICS,  NACA  0012  k  -  0  130  a 


(RADIANS) 


MOMENT  CHARACTERISTICS,  NACA  0012,  k  =  0.045,amav  =  30°,  U.=  41.7  fps. 


(SNVIQVH) 


MOMENT  CHARACTERISTICS, iNACA  0012,  k  =  0.079,a-  =  34°,  UU=  33.2  fps. 
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FIGURE  EM3 

COMPARISON  OF  LIFT  AND  MODEL  MOTION  TIME  HISTORY  -  EXAMPLE 
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FIGURE  ET-14 

DESCRIPTION  OF  "NOTCH"  DEPRESSION  IN  TIME  EVOLUTION  OF  LIFT  COEFFICIENT 


FIGURE  H-15 

VARIATION  OF  "NOTCH"  BEHAVIOR  WITH  DIMENSIONLESS  FREQUENCY 


V.  PARAMETRIC  COMPARISONS 


To  compare  airfoil  performance  over  the  broad  range  of  flow  and 
motion  conditions  examined  in  this  preliminary  study,  it  is  useful  to 
exploit  the  common  parameters  which  appear  in  the  test  matrix.  This 
section  of  the  report  will  focus  on  several  key  performance  "indicators" 
and  explore  their  variation  with  motion  amplitude,  dimensionless  pitch 
rate,  and  Reynolds  number. 

V.l  BEHAVIOR  OF  THE  SUCTION  PEAK 

The  influence  and  behavior  of  the  suspected  dynamic  stall  vortex 
can  be  inferred  by  examining  the  evolution  of  the  airfoil  upper  surface 
suction  peak.  One  description  of  this  behavior  is  provided  by  charting 
the  trajectory  of  this  minimum  pressure  location  as  a  function  of  the 
instantaneous  angle  of  attack.  These  results  are  shown  in  Figures  V-l 
through  V-3  for  three  different  values  of  the  dimensionless  frequency 
for  which  a  comparison  between  airfoils  was  possible. 

The  lowest  rate  case  displayed,  k  *  0.045,  is  also  the  minimum  in 
the  test  matrix  for  which  a  characteristic  vortex-induced  behavior  was 
observed.  Although  unsteady  effects  (including  time  lags  in  the  distribu¬ 
tions)  were  observed  at  much  lower  values  of  a  ,  an  identifiable  suction 
peak  which  "detached"  from  the  leading  edge  location  was  not  observed 
under  these  conditions.  As  Figure  V-l  shows,  the  suction  peak  location 
varies  linearly  with  the  incidence  angle  and,  therefore,  time,  over  a 
portion  of  the  curve  and  appears  to  be  independent  of  airfoil  type  in 
this  region.  For  the  specific  cases  described,  maximum  lift  coefficients 
were  observed  to  occur  in  the  23-28  degree  incidence  range  which  approxi¬ 
mately  coincides  with  the  onset  of  nonuniformities  in  the  suction  peak 


trajectories  in  the  three  curves.  In  both  NACA  0012  tests  displayed 
(different  Reynolds  numbers),  the  streamwise  rate  of  passage  of  the 
suction  peak  is  observed  to  decrease  as  the  maximum  lift  point  is  reached. 
In  contrast,  the  NACA  64iA012(13)  airfoil  exhibited  a  rather  unusual 
behavior  characterized  by  an  initial  decrease  in  suction  peak  longitudinal 
velocity  followed  by  a  rapid  increase  as  the  incidence  angle  continued 
to  increase. 

At  a  higher  dimensionless  pitch  rate  (k  =  0.09),  Figure  V-2  reveals 
that  the  initial  onset  of  suction  peak  motion  occurs  at  a  higher  incidence 
angle  when  compared  to  the  previous  results.  The  residence  of  the  suction 
peak  near  the  airfoil  leading  edge  is  reminiscent  of  steady  flow  behavior 
at  unstalled  incidence  angles.  Once  suction  peak  movement  in  the  down¬ 
stream  direction  begins,  it  proceeds  at  an  almost  constant  rate.  At 
this  point,  a  major  departure  from  the  behavior  observed  in  the  lower 
rate  case  occurs.  Over  a  range  of  angles,  the  location  of  the  suction 
peak  is  observed  to  remain  nearly  stationary  as  the  incidence  angle 
continues  to  increase.  This  result,  suggestive  of  a  "trapped"  vortex 
condition,  is  observed  for  both  airfoils,  although  the  conditions  for 
which  it  occurs  are  somewhat  different.  As  the  incidence  angle  increases 
further,  the  suction  peaks  once  again  begin  to  traverse  in  the  downstream 
direction  at  a  nearly  constant  rate  until  the  maximum  amplitude  of  motion 
is  reached.  It  is  useful  to  note  that  for  both  of  these  experiments, 
the  maximum  lift  coefficient  occurred  within  the  band  of  stationary 
suction  peak  behavior. 

At  the  highest  pitch  rate  shown  (Figure  V-3),  another  variation  in 
behavior  is  observed.  The  initial  departure  of  the  suction  peak  from 
the  leading  edge  region  occurs  at  yet  higher  incidence  angles  and  in  a 
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more  nonlinear  trajectory  than  for  the  lower  rate  cases.  In  both  examples 
shown  here,  the  suction  peak  traversal  speed  decreases  with  increasing 
incidence  angle  and,  therefore,  time.  Once  again,  a  stationary  suction 
peak  is  found  but  at  different  chordwise  locations  for  the  two  airfoils. 
This  location  is  near  the  30  per  cent  chord  point  for  the  NACA  0012 
airfoil  and  occurs  between  40  and  50  degrees  in  angle  of  attack.  For 
the  NACA  64]^A012(13)  airfoil  model,  the  suction  peak  remains  nearly 
stationary  at  or  near  the  20  per  cent  chord  point  between  32  and  44 
degrees  angle  of  attack  at  a  comparable  dimensionless  rate  of  motion. 

This  result  is  not  unexpected  given  the  predicted  influence  of  the 
pitching  axis  on  angle  at  which  stall  occurs.^  Subsequent  motion  of  the 
peak  downstream  is  observed  to  be  rapid  for  both  airfoils  up  to  the 
maximum  incidence  angle.  Comparisons  in  this  range  should  be  made 
cautiously,  however,  since  high  deceleration  rates  were  experienced  near 
the  maximum  amplitude  condition. 

At  this  point,  it  is  important  to  realize  that  observed  trends  and 
behavioral  variations  such  as  those  described  above  are  somewhat  affected 
by  the  spatial  resolution  of  the  measurements.  This  situation  is  a 
direct  result  of  the  limited  number  of  transducers  and  distribution  of 
pressure  measurement  locations.  For  example,  what  may  appear  as  an 
absolutely  stationary  phenomenon  (as  suggested  by  the  curves  in  the 
preceding  figure)  is  probably  dynamic  but  not  readily  resolved  due  to 
the  local  transducer  spacing  in  the  region  of  interest.  However,  the 
"relative"  stationarity  of  these  results  compared  to  other  portions  of 
the  cycle  is  discernable.  The  use  of  the  term  "stationary"  in  ensuing 
paragraphs  is  employed  in  this  relative  sense  recognizing  the  constraints 
imposed  by  the  limited  spatial  resolution  of  the  measurements. 


The  effects  of  maximum  motion  amplitude  (incidence  angle)  on  the 
suction  peak  evolution  at  fixed  dimensionless  pitch  rates  are  displayed 
in  Figures  V-4  and  V-5.  In  Figure  V-4,  this  variation  for  a  maximum 
incidence  angle  of  34  degrees  is  compared  with  the  very  large  amplitude 
case  described  previously.  A  mid-range  value  of  the  dimensionless  pitch 
rate  is  considered  here.  The  close  agreement  between  the  two  curves  at 
low  values  of  incidence  angle  is  testimony  to  the  repeatability  of  motion 
and  accuracy  of  the  measurements  in  these  experiments.  Departures  between 
the  two  trajectories  begin  to  occur  when  the  mrximum  angle  of  incidence 
is  reached  in  the  lower  amplitude  case.  However,  it  cannot  be  inferred 
that  vortex  entrapment  does  not  occur  in  this  case,  because  the  abscissa, 
a,  does  not  continue  to  vary  linearly  in  time  once  the  maximum  amplitude 
of  motion  is  reached.  Similar  behavior  is  documented  for  the  NACA 
64]A012(13)  airfoil  at  a  comparable  dimensionless  pitch  rate  as  shown  in 
Figure  V-5. 

In  addition  to  examining  the  trajectories  of  the  suction  peaks,  it 
would  be  useful  to  quantify  the  effectiveness  of  the  separation  vortex 
using  the  magnitude  of  the  suction  peak  as  it  traverses  the  surface. 

These  results  are  displayed  in  Figures  V-6  through  V-14.  In  these 
diagrams,  the  minimum  pressure  coefficient  is  plotted  against  the  instan¬ 
taneous  suction  peak  location.  Figure  V-6  illustrates  this  variation 
for  the  NACA  0012  airfoil  at  high  incidence  parameterized  by  different 
values  of  the  dimensionless  pitch  rate.  Note  that  while  the  magnitude 
of  the  suction  peak  increases  substantially  with  increased  pitch  rate, 
its  subsequent  decay  rate  is  also  commensurately  greater.  By  the  time 
the  suction  peak  has  reached  the  40  per  cent  chord  location,  the  pressure 
coefficient  has  increased  to  an  almost  constant  value  of  -3.0.  Subsequent 
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variations  of  this  parameter  as  the  pressure  peak  traverses  downstream 
follow  a  similar  trend  for  all  three  cases. 

The  insensitivity  of  the  separation  flowfield  to  airfoil  type  is 
perhaps  best  illustrated  in  Figures  V-7  through  V-9  in  which  a  comparison 
of  the  two  models  is  provided  for  three  different  dimensionless  pitch 
rates.  All  three  cases  involve  motion  to  very  large  incidence  angles 
(a  *  60  degrees).  While  some  differences  are  obvious  between  the  two 
airfoils  at  a  given  pitch  rate,  these  differences  are  small  compared  to 
the  variations  that  occur  as  a  function  of  pitch  rate. 

One  observation  apparent  from  Figure  V-7  is  that  the  suction  peak 
retains  its  "identity"  further  downstream  on  the  NACA  64  series  airfoil 
than  on  the  NACA  0012  airfoil. 

Figure  V-8  displays  comparable  results  for  a  higher  value  of  reduced 
frequency  again  showing  similar  trends  between  the  two  airfoils.  A 
noticeable  difference  from  the  lower  rate  case  occurs  for  both  airfoils 
near  the  30  per  cent  chord  point.  Coincident  with  the  stationary  vortex 
behavior  described  earlier,  the  suction  peak  amplitude  of  the  NACA 
64]A012(13)  airfoil  is  observed  to  increase  by  nearly  20  per  cent.  This 
suggests  a  strengthening  of  the  separation  vortex  during  that  period 
when  it  remains  nearly  stationary.  The  NACA  0012  airfoil  also  exhibits 
an  unusual  behavior  at  this  point  which  is  apparent  in  the  pressure 
coefficient  amplitude  time  history.  Specifically,  an  initial  decrease 
in  peak  amplitude  is  followed  by  a  small  increase  which  is,  in  turn, 
followed  by  another  decrease  as  the  suction  peak  continues  to  traverse 
downstream. 

In  Figure  V-9,  even  higher  pressure  coefficient  magnitudes  are 
observed,  but  the  increase  in  suction  peak  amplitude  at  the  stationary 


V-6 


vortex  location  is  not  evident.  It  is  apparent  from  these  figures  and 
those  described  above  that  the  dynamical  behavior  of  the  separated  flow 
region  changes  markedly  as  the  pitch  rate  is  varied.  It  further  suggests 
that  the  concept  of  identifying  the  local  surface  suction  peak  with  a 
single,  colocated  dynamic  stall  vortex  is  probably  an  oversimplified 
model  of  the  actual  physical  situation. 

Figures  V-10  through  V-12  provide  more  detailed  results  on  the 
evolution  of  suction  peak  magnitudes  with  location  for  the  NACA  64iA012(13) 
airfoil  model.  In  these  cases,  the  effects  of  maximum  incidence  angle 
on  flowfield  evolution  are  addressed.  For  the  lowest  rate  case  (k  = 

0.045)  shown  in  Figure  V-10,  the  three  sets  of  results  appear  very 
similar.  This  suggests  that  natural  separation  occurs  at  a  value  of 
lower  than  the  minimum  value  of  afflax  for  which  data  were  obtained.  The 
suction  peak  is  observed  to  reach  its  maximum  intensity  at  a  chordwise 
location  close  to  the  leading  edge  in  all  three  cases. 

A  similar  result  is  found  for  a  somewhat  higher  pitch  rate  as  shown 
in  Figure  V-ll.  Of  interest  here  is  the  fact  that  the  reintensification 
of  the  suction  peak  occurs  at  the  same  chordwise  location  and  is  of 
approximately  the  same  relative  magnitude  independent  of  the  maximum 
incidence  angle.  Even  after  the  peak  is  reached,  the  variation  in  the 
suction  peak  amplitude  between  the  two  curves  is  not  significant. 

An  even  stronger  similarity  in  behavior  is  apparent  in  Figure  V-12 
which  shows  the  same  result  for  the  highest  pitch  rate  case.  Again,  the 
three  sets  of  results  seem  to  point  at  the  dimensionless  pitch  rate  as 
being  a  singularly  important  parameter  which  dominates  the  physical 
development  of  these  flows. 
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V.  2  LIFT  AND  MOMENT  PEAKS 

An  assessment  of  unsteady  airfoil  performance  characteristics  is 

difficult  given  the  number  and  range  of  influential  performance  parameters 

and  the  additional  complication  of  arbitrary  airfoil  motion  time  history. 

In  this  section,  several  possible  measures  of  lift  enhancement  related 

to  separation  vortex  entrapment  are  proposed  and  described.  Resulting 

parametric  comparisons  using  these  variables  are  discussed  in  detail. 

Because  of  the  significance  of  moment  stall  effects,  similar  results  for 

peak  negative  moment  coefficient  will  also  be  discussed. 

The  simplest  measures  of  airfoil  performance  evident  from  force  and 

moment  time  histories  are  a  class  of  variables  which  can  be  termed  "peak 

parameters"  since  they  refer  to  the  maximum  and/or  minimum  conditions 

which  occur  during  the  motion  cycle.  These  include  but  are  not  limited 

to  the  angle  of  attack  for  which  the  lift  coefficient  is  a  maximum,  o^  , 

^max 

the  maximum  lift  coefficient,  Ci  ,  and  comparable  measures  of  the 

IDclX 

minimum  in  the  moment  coefficient,  namely  otr  and  C_  .  .  The  maximum 

^nin  “mm 

lift  coefficient  is  used  as  a  measure  of  the  severity  of  the  dynamic 

stall  and  the  influence  of  the  separation  vortex'  on  airfoil  lift.  The 

angle  of  attack  at  which  the  lift  coefficient  achieves  its  maximum  is 

hypothesized  to  be  associated  with  a  well-developed  separation  vortex, 

although  the  absence  of  detailed  flow  field  data  precludes  a  complete 

verification  of  this  assumption.  Similarly,  the  angle  otr  is 

^"min 

considered  representative  of  a  separation  process  which  has  matured, 
since  the  dynamic  stall  model  associates  this  event  with  the  passage  of 
the  stall  vortex  over  the  trailing  edge  region  of  the  airfoil.  The 
minimum  moment  coefficient,  C,^^,  can  be  viewed  as  a  measure  of  vortex 
organization  during  the  latter  stages  of  the  separation  process. 


The  variation  of  a0.  as  a  function  of  dimensionless  pitch  rate 

^max 

for  the  two  airfoils  is  depicted  in  Figure  V-13.  The  differences  in  the 

test  ranges  of  the  two  models  are  readily  apparent.  Results  have  been 

provided  for  several  values  of  the  maximum  amplitude  of  motion,  amax. 

Figure  V-13a  illustrates  that,  at  lower  pitch  rates,  “g.  increases  in 

1max 

linear  fashion  virtually  independent  »f  the  maximum  incidence  involved 

in  the  airfoil  incidence  time  history.  At  higher  pitch  rates  and  for 

lower  values  of  ct^ay,  the  angle  for  maximum  lift  coefficient  and  the 

maximum  amplitude  of  motion  are  nearly  coincident.  Even  for  extremely 

large  amplitude  motions,  the  value  of  a0.  appears  to  level  off  at 

xmax 

approximately  30  degrees  for  the  NACA  0012  airfoil. 

A  comparable  plot  for  the  NACA  64^A012(13)  airfoil  indicates  a 
similar  trend  near  a  value  of  k  ■  0.1,  as  shown  in  Figure  V-13b.  Proceeding 
to  higher  pitch  rates,  however,  revealed  a  resumption  in  the  increase  of 
the  maximum  lift  angle  with  pitch  rate.  This  behavior  is  somewhat 
anomolous  but  is  not  believed  to  be  uniquely  related  to  wind  tunnel 
characteristics  or  other  test  environment  idiosyncracies.  This  behavior 
is  recommended  as  a  subject  of  future  investigation.  A  comparison  of 
both  airfoils  in  the  region  of  measurement  overlap  shows  good  agreement 
between  the  two  sets  of  results. 

The  above  results  indirectly  suggest  that  the  onset  of  unsteady 
vortex  generation  might  be  affected  more  significantly  at  low  values  of 
the  dimensionless  pitch  rate  than  at  high  values.  That  is,  a  point  of 
diminishing  returns  is  soon  reached  for  lift  enhancement  as  one  proceeds 
to  higher  values  of  pitch  rate. 

The  maximum  lift  coefficient,  Ci  ,  plotted  as  a  function  of  the 
dimensionless  pitch  rate  is  depicted  for  the  two  airfoils  in  Figure 


V-14.  The  overall  shape  of  the  curves  is  remarkably  similar  to  that 
found  for  a  corresponding  incidence  angle  variation.  The  maximum  lift 
coefficient  increases  rapidly  and  monotonically  at  low  pitch  rates  with 
a  decrease  in  the  rate  of  rise  as  a  value  of  k  =  0.1  is  approached.  The 
magnitude  of  the  maximum  lift  coefficient  is  observed  to  be  dependent  on 
the  maximum  incidence  angle  of  airfoil  motion.  This  result  also  suggests 
that  the  rate  parameter  is  not  the  only  variable  which  can  influence 
airfoil  lift  performance. 

The  results  depicted  for  the  NACA  64]A012(13)  airfoil  in  Figure 
V-14b  confirm  that  the  growth  of  C].max  with  increased  pitch  rate  is  not 
as  large  at  high  values  of  k  as  it  is  for  lower  values.  Finally,  it  is 
instructive  to  compare  results  for  both  airfoils.  Despite  the  small 
range  of  overlap,  striking  similarities  in  gross  airfoil  performance  as 
measured  by  the  maximum  lift  coefficient  are  apparent. 

As  indicated  earlier,  the  strength  of  the  unsteady  separation  vortex 
near  the  end  of  the  useful  cycle  can  perhaps  best  be  characterized  by 
the  magnitude  of  the  minimum  pressure  peak,  Cmm£n-  The  variation  in 
this  parameter  as  a  function  of  the  dimensionless  pitch  rate  for  both 
airfoils  is  depicted  in  Figure  V-15.  From  the  results  of  Figure  V-15a, 
there  appears  to  be  a  significant  relationship  between  the  severity  of 
the  moment  stall  and  the  maximum  amplitude  of  motion,  amax  for  the  NACA 
0012  airfoil.  At  higher  pitch  rates,  more  severe  moment  stall  appears 
to  occur  at  yet  higher  values  of  incidence  angle.  In  cases  of  low 
maximum  incidence  angle,  those  for  which  separation  would  otherwise 
occur  at  values  higher  than  the  actual  maximum  incidence  of  motion,  the 
minimum  peak  pressure  coefficient  is  actually  observed  to  increase 
(magnitude  of  minimum  pressure  spike  decreases)  with  increasing  pitch 
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A  comparable  result  for  the  NACA  64^A012(13)  airfoil  shown  in  Figure 
V-15b  appears  anomolous,  especially  for  the  largest  incidence  angle 
depicted.  However,  the  trend  to  higher  suction  peaks  with  increasing 
motion  amplitude  is  again  evident  except  for  the  very  highest  pitch 
rates  measured.  A  comparison  of  the  two  airfoils  in  the  overlap  region 
again  shows  good  agreement. 

The  anomolous  behavior  in  minimum  pressure  coefficient  at  the  highest 
reduced  frequencies  has  not  been  adequately  explained  by  any  companion 
measurements,  and  remains  a  topic  for  future  investigation.  However,  an 
interesting  observation  can  be  made  by  comparing  the  variation  of  maximum 
lift  coefficient  with  that  of  the  minimum  moment  coefficient;  While 
both  parameters  initially  increase  with  pitch  rate  (or  frequency),  the 
maximum  lift  coefficient  appears  to  reach  a  plateau  at  a  lower  value  of 
the  dimensionless  pitch  rate  than  the  corresponding  minimum  moment  coeffi¬ 
cient.  The  lack  of  correlation  of  these  results  is  also  recommended  as 
a  subject  of  further  study. 

V. 3  AIRFOIL  PERFORMANCE  COMPARISONS  -  COMPLETE  MOTION  CYCLE 

Although  the  peak  values  of  the  variables  described  in  the  preceding 
section  provide  a  gross  comparison  between  the  various  test  cases,  they 
do  not  fully  account  for  the  effects  of  the  model  motion  time  history  on 
overall  airfoil  performance.  It  is  useful,  therefore,  to  devise  other 
parameters  which  can  more  completely  characterize  the  integral  performance 
effects  over  the  motion  cycle.  It  would  also  be  useful  to  devise  a 
parametric  scheme  which  would  allow  for  comparison  of  airfoil  performance 
between  cases  where  the  motion  time  history  is  fundamentally  different. 

In  examining  the  various  force  and  moment  coefficient  time  histories, 
the  most  obvious  parameter  to  be  included  in  any  such  characterization 
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is  Che  time  itself.  In  addition,  the  time  associated  with  any  given 
event  can  be  nondimens ionalized  by  a  characteristic  time  which  can  either 
be  associated  with  the  model  performance  or  one  related  to  the  reference 
flow  conditions.  If  a  nondimens ionalizat ion  independent  of  the  model 
motion  time  history  is  desired,  a  reference  time  scale  based  on  external 
flow  conditions  must  be  employed.  In  this  case,  t  -  c/U^,  is  selected, 
where  t  is  the  time  it  takes  for  a  freestream  fluid  element  to  travel  a 
distance  equal  to  the  length  of  the  airfoil  chord. 

Another  important  question  is  concerned  with  the  definition  of  an 
appropriate  time  measurement  which  characterizes  the  effectiveness  of  a 
given  motion  cycle.  There  appears  to  be  no  single  best  answer  to  this 
question  since  the  resulting  measurement  depends  to  some  degree  on  the 
element  of  performance  to  be  assessed.  Since  this  report  deals  primarily 
with  the  optimization  of  airfoil  lift  characteristics,  the  time  that 
defines  the  persistence  of  lift  coefficient  above  the  static  stall  value, 
t+  (as  shown  in  Figure  IV-13),  will  be  employed  here.  Referenced  to  the 
freestream  variable,  t,  this  temporal  parameter  is  hereafter  referred  to 


as  the  duration  coefficient,  C^,  where 

c  =  t+ 

The  variation  of  this  parameter  with  the  dimensionless  pitch  rate 
for  the  two  airfoils  is  shown  in  Figure  V-16.  In  part  (a)  of  this  figure, 
results  for  the  NACA  0012  airfoil  model  are  depicted  for  several  different 
pitch  amplitudes.  The  symbols  are  used  to  characterize  results  which 
are  within  +  3  degrees  of  a  nominal  value.  These  results  can  then  facili¬ 
tate  a  comparison  of  the  two  airfoils  under  comparable  conditions.  The 
increase  in  duration  coefficient  occurs  rapidly  and  monotonically  for 
low  values  of  the  dimensionless  pitch  rate  (k  <  0.015).  A  single  anomolous 
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data  point  occurred  at  a  k-value  of  0.0182  where  Ct  was  measured  to  be 
17.2.  The  potential  error  in  this  measurement  is  estimated  to  be  high 
due  to  the  shallow  negative  slope  of  the  C\  -  time  history  curve  near 
the  static  stall  point.  The  remaining  data  indicate  that,  in  the  mid¬ 
range  region,  the  rate  of  increase  slows  until  a  nearly  constant  value 
is  achieved  over  the  highest  set  of  pitch-rates  examined.  The  effect  of 
maximum  pitch  amplitude  on  this  characteristic  cannot  be  determined  with 
any  degree  of  accuracy  due  to  the  sparcity  of  available  data.  However, 
the  evidence  suggests  that  the  maximum  value  of  the  incidence  angle  for 
constant  rate  and  hold  motion  does  play  a  role  in  the  determination  of 
this  coefficient. 

Results  for  the  NACA  64iA012(13)  airfoil  are  shown  in  Figure  V-l6b. 
Although  these  data  cover  a  wider  range  of  reduced  pitch  rates,  insuffi¬ 
cient  results  are  available  at  the  lowest  values  to  assess  trends  in 
that  region.  However,  the  curve  representing  the  highest  amplitude 
motion  (amax  =  60°)  does  provide  yet  another  interesting  insight.  In 
contrast  with  previous  trends,  the  duration  coefficient  is  observed  to 
decrease  with  continued  increases  in  the  pitch  rate.  This  is  not  the 
case  at  the  lower  value  of  the  maximum  incidence  angle  plotted  in  the 
figure.  More  significantly,  measured  values  of  the  duration  coefficient 
for  the  40  degree  case  are  higher  than  their  counterparts  for  the  higher 
maximum  incidence  angle  at  comparable  values  of  the  reduced  pitch  rate. 

The  single  data  point  plotted  for  the  lowest  rate  suggests  that  a  decrease 
of  an  additional  10  degrees  in  the  value  of  amax  results  in  a  much  decreased 
value  of  the  duration  coefficient.  These  combined  results  suggest  that 
airfoil  performance  improvements  do  not  necessarily  correlate  directly 
with  increases  in  the  pitch  amplitude  of  motion,  but  that  an  optimum 


value  of  %ax  exists  and  probably  depends  on  the  pitch  rate  and  other 
flow  conditions.  This  inference  will  be  examined  more  completely  later 


in  this  report. 

While  the  duration  coefficient  provides  a  useful  relative  measure 
of  performance  comparison,  it  is  susceptible  to  large  inherent  errors 
due  to  the  lack  of  repeatability  of  deterministic  flow  processes  (stall 
vortex  coherence)  near  the  end  of  the  motion  cycle.  Under  some  conditions, 
the  statistical  scatter  in  the  measurement  of  t+  can  be  unacceptably 
high.  Moreover,  the  duration  coefficient  does  not  account  for  the  magnitude 
of  lift  overshoot  or  its  variation  in  time.  A  more  comprehensive  indicator 
must,  therefore,  incorporate  these  effects  in  some  integral  fashion.  An 
obvious  candidate  that  satisfies  these  criteria  and  which  is  also  independent 
of  the  type  of  motion  is  the  area  under  the  force-time  curve  which  denotes 
an  effective  additional  impulse,  II-  Since  lift  optimization  is  a  principal 
objective  of  this  investigation,  the  parameter  is  here  defined  for  values 
of  lift  greater  than  the  maximum  steady  flow  value.  In  dimensional 
form, 

IL  =  J  4L-.it 

o 

for  all 

“■  s  i  o 

A  non-dimensional  version  of  this  relationship  yielding  a  lift 
impulse  coefficient  can  be  written,  in  general,  as 
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for  all 
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Ci  >  0 
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where  tn  is  an  appropriate  reference  time  scale. 

Two  non-dimensionalization  schemes  (reference  time  scales)  are 
physically  significant.  If  tn  is  taken  as  the  total  time  for  which  the 
lift  exceeds  the  maximum  steady  flow  value,  then, 


“n 


The  impulse  relation  then  becomes 

ACj 


/  W  • 


(9) 


for  all 

4Ci 5  °i  -  °w8  ±  0 

In  this  case,  the  result  is  equal  to  the  average  value  of  the  excess 
lift  increment. 

The  behavior  of  this  average  excess  lift  parameter  for  the  NACA 
0012  airfoil  is  documented  in  Figure  V-17.  The  upper  (a)  portion  in 
this  figure  illustrates  the  variation  with  dimensionless  pitch  rate  for 
the  highest  motion  amplitudes  examined  (amax  =  60°).  The  effect  of 
motion  amplitude  is  documented  separately  in  Figure  V-17(b). 

The  average  excess  lift  coefficient  is  observed  to  vary  with  pitch 
rate  in  a  manner  not  unlike  that  documented  for  the  section  pressure 
peak  described  earlier.  A  rapid  increase  with  pitch  rate  at  low  values 
of  that  parameter  are  followed  by  a  gradual  flattening  of  the  curve 
until  a  nearly  constant  level  is  achieved  at  the  highest  levels.  A 
least  squares  curve  is  fitted  through  the  available  data.  Reynolds 
number  independence  ovar  the  ran  'e  of  study  is  evident  by  the  similarity 


in  results  achieved  ai 


&l\  of  k  ■  0.045. 


The  variation  of  maximum  incidence  angle  with  the  average  incre¬ 
mental  lift  coefficient  (Figure  V-17b)  shows  that,  at  least  for  higher 
rates  of  motion,  the  optimum  incidence  angle  lies  above  the  static  stall 
value,  but  not  at  the  highest  angles  examined.  It  appears  that  optimum 
lift  performance  might  be  achieved  when  the  maximum  incidence  angle  of 
the  motion  is  comparable  to  the  expected  separation  angle  associated 
with  the  pitch  rate  employed. 

A  similar  result  for  the  NACA  64]A012(13)  airfoil  is  provided  in 
Figure  V-18.  Both  pitch  rate  and  maximum  amplitude  effects  have  been 
condensed  to  a  single  set  of  curves  in  this  case.  Several  overlap  points 
from  the  NACA  0012  data  have  been  plotted  to  compare  the  two  airfoils. 

The  insensitivity  of  airfoil  type  under  these  motion  extremes  is  again 
evident.  Note  that  the  excess  lift  is  nearly  constant  with  a  slight 
downward  trend  over  this  range  of  high  pitch  rates.  In  addition,  it  is 
very  apparent  that  the  lower  value  of  amax  results  in  consistently  higher 
performance  in  these  instances.  It  is  also  evident  that  a  larger  number 
n*  test  cases  is  desirable  if  the  excess  lift  is  to  be  maximized. 

Recalling  the  non-dimensional  impulse  relation,  eqn.  (8),  it  is 
useful  to  consider  a  different  reference  time  scale  based  on  the  free- 
stream  flow  conditions,  i.e., 


Then,  we  can  write 


f  *C1 


d<T> 


(10) 


/s 
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The  variable,  Ic,  will  be  referred  to  as  the  "dimensionless  impulse." 

This  parameter  is  a  versatile  measure  of  performance  since  it  incorporates 
both  lift  magnitude  and  duration.  In  addition,  it  is  capable  of  providing 
equitable  comparisons  of  highly  dissimilar  motion  time  histories  since 
the  freestream  time  scale  is  accounted  for  implicitely. 

The  variation  of  Ic  with  dimensionless  pitch  rate  for  the  NACA  0012 
model  is  shown  in  Figure  V-19a.  Results  are  depicted  for  several  common 
and  approximately  equal  (  +  3°)  values  of  the  maximum  pitch  incidence 
angle.  The  characteristic  shape  of  the  maximum  incidence  (amax  =  60°) 
characteristic  is  familiar  by  now.  A  rapid  increase  in  the  impulse 
parameter  at  low  pitch  rates  is  followed  by  a  more  gradual  increase  at 
the  higher  rates  of  motion.  The  tendency  of  the  curve  to  "flatten"  and 
achieve  a  local  maximum  is  not  apparent  in  the  behavior  of  this  parameter. 

The  limited  results  available  for  other  values  of  maximum  pitch 
amplitude  preclude  a  detailed  discussion  of  their  behavior  in  the  context 
of  this  figure.  However,  a  limited  comparison  of  motion  amplitude  can 
be  effected  through  the  crossplots  provided  in  Figure  V-19b.  From  these 
curves,  one  observes  that  the  incidence  angle  selected  for  maximum  perfor¬ 
mance  is  very  much  dependent  on  the  pitch  rate  employed. 

A  comparable  set  of  results  for  the  NACA  64iA012(13)  airfoil  is 
provided  in  Figure  V-20.  For  this  wider  range  of  reduced  pitch  rates, 
several  contradictory  trends  are  observed.  Tests  at  high  incidence 
angle  (<*niax  *  60°)  resulted  in  a  reduction  in  the  dimensionless  impulse 
with  increased  values  of  the  reduced  pitch  rate  (on  the  average)  as 
shown  in  Figure  V-20a.  A  nearly  opposite  trend  is  observed  when  the 


maximum  incidence  of  motion  approximates  the  angle  at  which  leading  edge 
separation  would  otherwise  occur  had  the  motion  continued  without  inter¬ 
ruption. 

The  significance  of  maximum  incidence  amplitude  is  even  more  apparent 
in  Figure  V-20b.  For  all  pitch  rates  plotted,  the  optimum  airfoil  perfor¬ 
mance  occurred  at  an  incidence  angle  of  approximately  40  degrees.  This 
coincides  with  values  close  to  the  observed  angle  for  leading  edge  separa¬ 
tion  (i.e.,  suction  peak  detachment)  in  this  pitch  rate  range. 

A  comparison  of  the  results  in  the  region  of  measurement  overlap 
shows  reasonable  agreement  at  a  reduced  pitch  rate  near  k  =  0.10.  How¬ 
ever,  there  is  a  substantial  disagreement  in  the  results  at  lower  pitch 
rates  (k  =  0.05).  The  cource  of  this  disagreement  warrants  further 
investigation. 

V.4  CORRELATION  OF  SUCTION  PEAK  BEHAVIOR  WITH  AIRFOIL  PERFORMANCE 

The  influence  of  the  motion  and  magnitude  of  the  suction  peak  (and, 
therefore,  of  the  unsteady  vortex)  on  airfoil  lift  performance  can  be 
assessed  through  several  parametric  cross  comparisons.  One  useful  measure 
is  provided  by  relating  the  suction  peak  trajectory  to  corresponding 
values  of  instantaneous  lift  coefficient.  These  results  are  provided 
for  large  amplitude  motion  cases  in  Figures  V-21  through  V-24. 

The  effect  of  dimensionless  pitch  rate  on  performance  of  the  NACA 
0012  airfoil  is  evident  in  the  plots  of  Figure  V-21.  For  the  lowest 
value  of  k,  the  maximum  lift  coefficient  is  observed  to  occur  when  the 
suction  peak  resides  at  a  chordwise  location  of  approximately  40  per 
cent.  For  the  intermediate  pitch  rate  case,  the  location  of  the  suction 
peak  for  maximum  lift  appears  to  move  forward  to  nearly  the  30  per  cent 
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chord  point.  The  unusual  reenergization  behavior  discussed  previously 
is  also  apparent  in  this  middle  curve.  Its  appearance  in  this  figure 
confirms  that  the  increase  in  suction  peak  intensity  during  the  stationary 
period  coincides  with  a  net  increase  in  airfoil  lift.  The  net  increase 
in  lift  during  this  static  phase  is  nearly  12  per  cent  in  the  present 
case. 

A  comparison  of  the  two  airfoils  is  provided  in  Figures  V-22  through 
V-24.  The  lowest  pitch  rate  case  of  Figure  V-22  illustrates  comparable 
behavior  until  the  suction  peak  reaches  the  30  per  cent  chord  location. 

At  this  point,  the  lift  of  the  NACA  64]A012(13)  is  observed  to  increase 
substantially  as  the  suction  peak  remains  near  stationary. 

Figure  V-23  describes  the  situation  for  a  higher  reduced  pitch  rate 
and  one  in  which  a  stationary  suction  peak  was  observed  for  both  airfoils. 
The  increase  in  peak  suction  coincides  with  an  increase  in  the  lift 
coefficient.  However,  the  percentage  change  in  the  lift  force  is  some¬ 
what  smaller  than  the  net  change  in  the  local  suction  peak  amplitude. 

A  similar  comparison  provided  in  Figure  V-24  for  the  highest  pitch 
rate  case  indicates  that,  although  the  shapes  of  the  two  curves  are 
similar,  the  characteristic  trajectory  is  slightly  shifted  between  the 
two  airfoils.  The  maximum  lift  occurs  earlier  in  the  evolution  of  the 
suction  peak  (separation  vortex)  for  the  64-series  airfoil  than  for  the 
NACA  0012.  These  discrepancies  may  be  attributable  to  the  difference  in 
pivot  (rotation)  location  between  the  two  airfoils  and  not  necessarily 
to  the  variations  in  airfoil  shape. 

The  influence  of  pitch  rate  on  the  suction  peak  location  for  which 
the  lift  is  maximum  is  most  easily  observed  in  a  consolidated  version  of 
the  previous  figures  which  is  provided  in  Figure  V-25.  It  is  apparent 


that  the  primary  effect  of  increasing  the  pitch  rate  for  any  given  set 
of  flow  conditions  is  to  cause  the  location  of  the  separation  vortex  for 
maximum  lift  to  move  nearer  to  the  leading  edge. 

Finally,  one  can  examine  the  variation  in  lift  coefficient  with  the 
magnitude  of  the  suction  peak  as  it  evolves  along  the  surface.  Curves 
illustrating  this  behavior  for  the  same  cases  described  in  the  preceding 
paragraphs  are  provided  in  Figures  V-26  through  V-30. 

These  results  for  the  NACA  0012  airfoil  are  overlaid  for  a  range  of 
reduced  pitch  rates  in  Figure  V-26.  The  correspondence  between  increasing 
lift  coefficient  and  increasing  amplitude  of  the  suction  peak  is  evident. 
The  hysteresis  behavior  in  these  time  histories  is  also  apparent.  Figure 
V-27  shows  a  comparison  between  the  NACA  0012  airfoil  and  the  NACA 
64]A012(13)  airfoil  at  nearly  identical  pitch  rates  and  very  high  values 
of  the  maximum  incidence.  As  has  been  the  case  in  most  cases  discussed 
previously,  the  behavioral  differences  between  the  two  airfoils  are 
small  compared  to  those  related  to  the  dimensionless  pitch  rate  variation. 
This  is  equally  apparent  in  Figures  V-28  and  V-29.  A  comparison  of 
these  high  amplitude  motion  cases  reveals  that  increases  in  the  pitch 
rate  result  in  correlatable  increases  in  both  the  suction  peak  magnitude 
and  the  overall  lift  coefficient.  That  is,  the  lifting  capability  of 
the  airfoil  under  these  large  amplitude  excursions  is  inextricably  linked 
to  the  strength  of  the  developing  unsteady  separation  vortex  manifested 
through  the  magnitude  of  the  suction  peak. 

The  influence  of  maximum  incidence  angle  at  nearly  comparable  values 
of  the  reduced  pitch  rate  (a  12  per  cent  variation,  in  this  case)  is 
displayed  in  Figure  V-30  for  the  NACA  0012  airfoil.  While  minor  discrep¬ 
ancies  appear  to  occur  following  the  onset  of  separation,  the  evolutionary 


V-20 


behavior  of  the  two  curves  is  very  similar  when  compared  to  similar 
curves  for  different  values  of  k. 

An  additional  correlation  of  these  results  is  possible  from  the 
condensed  format  provided  in  Figure  V-31.  The  nearly  linear  increase 
with  pitch  rate  confirms  the  relationship  between  suction  peak  and  lift 
coefficient  and  reinforces  the  notion  that  separation  vortex  strength 
plays  an  important  role  in  determining  airfoil  performance  character¬ 
istics  after  separation  occurs. 
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NACA  0012,  a  =  60° 


VI.  SUMMARY  AND  CONCLUSIONS 


Results  of  the  present  investigation  involving  those  measurements 
which  focus  on  the  effects  of  constant  pitch  rate  airfoil  motions  can  be 
summarized  as  follows: 

(a)  The  generation  of  airfoil  constant  pitch  rate  motions  over  a 
wide  range  of  parametric  flow  conditions  has  been  successfully  demon¬ 
strated  using  the  FJSRL  two-degree-of- freedom  unsteady  motion  apparatus. 
Cycle-to-cycle  motion  repeatability  was  observed  to  be  excellent. 

(b)  For  both  airfoils  investigated,  increases  in  pitch  rate  were 
found  to  further  delay  the  onset  of  unsteady  leading  edge  separation  as 
inferred  from  upper  airfoil  surface  pressure  distributions;  However, 
the  delay  was  found  to  increase  less  rapidly  at  progressively  higher 
values  of  the  dimensionless  pitch  rate. 

(c)  At  incidence  angles  well  in  excess  of  the  static  stall  value, 
the  appearance  of  a  steep  suction  spike  on  the  airfoil  upper  surface  was 
followed  by  its  subsequent  reduction  in  intensity  and  spatial  diffusion 
as  it  traversed  downstream.  This  phenomenological  behavior,  apparent  to 
varying  degrees  in  all  of  the  experiments,  is  consistent  with  that 
observed  by  other  investigators  principally  in  harmonically-induced 
dynamic  stall  experiments.  The  evolution  of  this  peak  is  believed  to  be 
controlled  by  the  behavior  of  the  unsteady  leading  edge  separation  vortex. 

(d)  Maximum  lift  coefficients  in  excess  of  three  times  the  steady 
flow  (no  motion)  static  lift  coefficient  were  observed  in  some  cases  at 
high  pitch  rates. 

(e)  A  dimensionless  lift  impulse  function  was  defined  and  proposed 
as  a  suitable  measure  for  comparing  the  effectiveness  of  various  motion 
types  on  sustained  lift  optimization. 
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(f)  It  was  observed  that  neither  the  highest  motion  rates  nor  the 
largest  pitch  amplitudes  resulted  in  the  maximum  values  of  the  pitch 
impulse  parameter.  For  rigid  airfoil  rotation  at  constant  rate,  the 
maximization  of  the  unsteady  lift  appears  to  occur  at  intermediate  values 
of  the  pitch  rate.  For  a  given  rate  (<*),  the  optimum  impulse  occurs  when 
the  maximum  incidence  angle  closely  approximates  the  angle  at  which 
leading  edge  separation  would  have  been  expected  had  the  pitch-up  motion 


continued  unaltered. 
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APPENDIX 


Table  III-l 


Model  Motion  Parameters  -  NACA  0012 


Pitch  Motion  Translational  Motion 


Test 

//  Uoo(fps) 

TyPe'1' 

k  =  ctc/2Uco* 

a<°> 

max 

Type+  Rate/Freq 

Notes 

1 

41.6 

A 

0.001 

31 

None 

Constant  pitch 

2 

41.2 

A 

0.001 

22 

rate  motion 

3 

41.2 

A 

0.005 

57 

only 

4 

41.3 

A 

0.005 

41 

5 

41.9 

A 

0.005 

22 

6 

41.9 

A 

0.009 

61 

7 

41.4 

A 

0.009 

41 

8 

41.4 

A 

0.010 

21 

9 

41.7 

A 

0.018 

60 

10 

33.1 

A 

0.047 

60 

11 

41.2 

A 

0.046 

60 

12 

67.9 

A 

0.044 

58 

13 

41.7 

A 

0.045 

30 

14 

40.9 

A 

0.073 

57 

15 

33.1 

A 

0.089 

56 

16 

33.2 

A 

0.079 

34 

17 

33.2 

A 

0.081 

25 

18 

33.2 

A 

0.080 

15 

19 

33.2 

A 

0.108 

55 

' 

r 

20 

33.9 

A 

0.130 

55 

1 

21 

41.4 

B 

0.042 

27 

None 

Constant  pitch 

22 

33.2 

B 

0.076 

23 

w/abrupt  return 

23 

33.6 

B 

0.101 

22 

1 

’ 

to  initial 

incidence  angle 

24 

41.5 

C 

a  *  .45  rad/s2 

56 

None 

Constant  pitch 

25 

41.2 

C 

<i  *  16.1  rad/s2 

62 

acceleration 

26 

41.9 

C 

a  *  66, 4  rad/s2 

61 

only 

27 

41.0 

C 

8  =  152.0  rad/s2 

60 

J 

28 

41.3 

C 

a  »  200.2  rad/s2 

61 

1 

t 

29 

42.6 

A 

k-  0.058 

27 

F 

x  *  19.7  ft/s2 

Combined  motion 

w/delayed  pitch 

up 

+see 

4.  . 

Figure  III-3 

*dlmens ionless  parameter,  k,  based  on  pitch  rate  unless  otherwise  noted 
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Table  I I I- 2 


Model  Motion  Parameters  -  NACA  64jA012(13) 


Pitch  Motion  Translational  Motion 


Test  # 

U»(fps) 

Type* 

k-dc/2U» 

Type"*" 

Rate/Freq 

Notes 

30 

32.1 

A 

0.045 

37 

None 

Constant  pitch  rate 

31 

32.1 

A 

0.045 

58 

motion  only 

32 

32.1 

A 

0.084 

37 

k  “OC/2U,, 

33 

32.1 

A 

0.094 

59 

34 

32.1 

A 

0.131 

38 

35 

32.5 

A 

0.136 

59 

36 

32.5 

A 

0.171 

39 

37 

32.8 

A 

0.173 

60 

38 

32.1 

A 

0.211 

59 

> 

39 

32.1 

B 

0.043 

29 

’ 

40 

32.1 

None 

12 

E 

x  *  -5.0  fps 

Airfoil 

at  fixed 

41 

32.1 

None 

20 

E 

x  ”  -5.0  fps 

incidence 

42 

32.5 

A 

0.137 

19 

E 

x  ”  -4.8  fps 

Combined  motion  w/ 

43 

32.1 

A 

0.143 

38 

E 

x  * -5.1  fps 

delayed  start  in  pitch 

44 

32.1 

A 

0.143 

20 

E 

x  ■  4.9  fps 

1 

45 

32.2 

A 

0.145 

39 

E 

x  *=  4.8  f ps 

1 

k'  -  0JC/2TJa, 

a£(°) 

46 

32.1 

D 

0.133 

5  None 

“max  *30°,  all  harmonic 

47 

32.1 

D 

0.134 

10 

cases 

48 

32.0 

D 

0.135 

16 

49 

32.0 

D 

0.135 

21 

50 

32.1 

D 

0.135 

26 

51 

32.1 

D 

0.244 

10 

52 

32.1 

D 

0. 246 

26 

• 

53 

32.1 

D 

0.373 

10 

54 

32.1 

D 

0.375 

26 

+see  Pigure  II1-3 
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ERRATA 


The  following  corrections  should  be  annotated  with  reference  to 
the  figures  and  illustrations: 


Page 

Figure  Number 

Correction 

11-21 

II-7 

An  additional  pressure  transducer 
should  appear  at  the  airfoil  leading 
edge,  x/c  =  0.0;  note  with  "x" . 

111-8,9 

111-1,2 

Axes  labels  missing;  Abscissa  label 
is  "x/c";  Ordinate  label  is  "Cp" . 

III-12  - 
III-30 

III-5  -  III-10 

Axes  labels  missing;  Abscissa  label 
is  "x/c" ;  Ordinate  label  is  "Cp" . 

III-16 

III-6 

Figure  (H)  label  incorrect;  should  be 
labeled  "34.4°". 

IV- 6 

IV- 2 

Change  label  from  "Steady  Flow  Lift 
Coefficients,  NACA  0012  Airfoil"  to 
"Steady  Flow  Lift  Coefficients,  NACA 
64^012(13)  Airfoil". 

IV-17 

IV-13 

Change  location  label  on  ordinate 

from  "C.  "  to  "C»  ". 

•max  *maxs 

IV- 18 

IV- 14 

Change  location  label  on  abscissa 

from  "t-.ax-C."  to  "tr. 

*  H  max 

IV-19 

IV- 15 

4  labels  in  legend  should  be  changed 
from  " . . . ,  max  ="  to  max  ="  . 

V-21 

V-l 

Ordinate  label  missing;  label  should 

read  "(x/c)p  . 

Fmin 

V-33 

V-13 (a) 

Label  "(A)"  should  read  "(a)";  ordinate 

label  "vcpmx"  should  read  "crCj 

a 


V-37,38  V-15 


Ordinate  label  "CM._ "  should  read 

Up  ••  WMN 

C«mn  • 


V-39 

V-16 (a) 

Zn  legend,  3rd  line  from  top  should 
note  "$  =  60°,  U*  =  68  fps". 

V— 41 

v-17 (a) 

Ordinate  label  "  AC , "  should  read 
”AC| " ;  in  legend,  6th  label  from  top 
should  be  changed  from  "O"  to  "®" . 

V-42 

V-17 (b) 

Figure  title  should  read  "Variation 
with  tfmax" * 

V-44 

V-19 (a) 

In  legend,  6th  symbol  from  top  should 
be  changed  from  "O"  to  "®". 

V-48  - 
V-51 

V-21  -  V- 24 

Abscissa  label  should  read  "(x/c)n 

^mm 

V-50 

V-2  3 

In  legend,  change  "NACA  64^012(13)1..." 
to  "NACA  64^012(13),  ...". 

V-51 

V-24 

In  figure  title,  change  "...  fps, 
max  “  to  "...  fps,  «max  =  60°-. 

V-58 

V-31 

In  figure  title,  change  "...  OR 
MAXIMUM..."  to  "...  AT  MAXIMUM. . .". 

